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ABSTRACT 



The flowfield downstream of a vertically-launched surface-to-air missile 
model at an angle of attack of 50° and a Reynolds number of 1.1 x 10 5 was 
investigated in a wind tunnel of the Naval Postgraduate School. The goal of this 
thesis is to experimentally validate the pressure measurement system for flowfield 
variables with elevated levels of turbulence; to determine the location and intensity 
of the asymmetric vortices in the wake of the VLSAM model at a raised level of 
freestream turbulence; and to display the asymmetric vortices by velocity mapping 
and pressure contours. The purpose is to correlate the results with the force 
measurements of Rabang to provide a greater understanding of the vortex 
flowfield. The body-only configuration was tested. Two flowfield conditions were 
treated: the nominal ambient wind tunnel condition, and a condition with grid- 
generated turbulence of 3.8% turbulence intensity and a dissipation length scale of 
1.7 inches. The following conclusions were reached: 1) The relative strengths of 
the asymmetric vortices can be noted by the sharp spike shape in the ambient 
condition; this condition becomes diffused and becomes fatter in the turbulent 
condition; 2) The right side vortex has greater strength than the left side one as 
seen by the diffusion in the total pressure coefficient and static pressure coefficient 
contours with and without a turbulent condition; 3) an increase in turbulence 
intensity tends to reduce the strength of the asymmetric nose-generated vortices; 
also pushes the two asymmetric vortices closer together; 4) and crossflow 
velocities were examined and were found to indicate the behavior denoted by the 
pressure contours. 
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NOMENCLATURE 



a: angle of attack 

Oav-' the AOA at which steady asymmetric vortices are formed 
AOA: angle of attack 

ots V : the AOA at which steady symmetric vortices are formed 

otuv- the AOA at which unsteady vortices are formed 

Cp S : Static pressure coefficient 

Cpt-' Total pressure coefficient 
C y : side force coefficient 

d: base diameter of the missile body 

e: blockage correction 

K: wind tunnel calibration factor 

Lj: missile length scale 

L^: missile diameter 

Le: turbulent scale 

l n /d: missile's nose fineness ratio 

l n : nose length 

M: mach number 

N: normal force 

Pf. indicated total pressure 

P2: indicated static pressure 

P3: indicated static pressure 

P4: pitch angle pressure 

P5: pitch angle pressure 
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P atm : ambient pressure 

P s : static pressure 

P t : total pressure 

Ptc: total pressure coefficient of pressure probe calibration 

0: pitch angle 

Q: free stream dynamic pressure 

q: local dynamic pressure 

0a: nose semi-vortex angle 

0c: pitch angle coefficient 

q M : measured dynamic pressure 

0r: roll angle 

R: air constant - 1718 

p: air density 

Re: Reynolds number C = (pUd/p.) 

T : temperature 

t: time 

T a v: average temperature 

T 0 : ambient temperature 

T u : turbulence intensity - (u/U) 

U: longitudinal mean velocity 

u: longitudinal velocity function 

U m : Test section velocity 

v: lateral velocity function 

V c : velocity coefficient of pressure probe calibration 

V v : vertical velocity 
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Y : side force 

ji: absolute viscosity 
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I. INTRODUCTION 



A. BACKGROUND 

In the past there has been a need for trainable launchers, both in azimuth and in 
elevation, so that the anti-aircraft missiles could be fired into gathering and 
guidance beams which direct them toward their targets. By way of contrast, a 
missile fired from a vertical launcher has to guide itself into the appropriate 
direction and heading. [Ref.l] 

Vertical launchers are part of a general evolution in anti-aircraft missile 
systems towards a higher rate of target engagement, which has three 
complementary aspects. The first aspect is the number of targets arriving 
simultaneously or nearly simultaneously. The ship must attempt to overcome the 
usual limitation of one engagement per fire control radar. Having overcome this 
limitation, the ship must keep enough defensive missiles in the air at any one time to 
fully exploit fire control capacity. The US AEGIS class ships are equipped with the 
US Mark 26 trainable launcher which fires about twice as fast as its predecessor, the 
Mark 13. The new Martin Mark 41 vertical launcher fires about five times faster 
than the US Mark 13 trainable launcher, i.e., at about one missile per second. 
[Ref. 1] 

The second aspect is the number of targets which appear over an extended 
length of time. Current point defense systems, such as SEA SPARROW and SEA 
WOLF, generally employ box launchers containing six or eight rounds. Once they 
have been fired off, the launcher must be reloaded by hand, a relatively laborious 
process. To provide automatic reloading would require considerable below deck 



1 



space, which is in short supply on most ships. A ship with a vertical launcher and 
with each missile stored in and launched from its own canister, offers a significant 
improvement in magazine handling within the usual constraints of above-deck 
space. 

The third aspect is the ability to carry out sustained operation in the face of air 
opposition. Missiles are voracious consumers of magazine space, hence the number 
of missiles per ship is always very limited. A Magazine/Launch system, then, has to 
be adapted to rapid replenishment if ships are to operate for a protracted period. 
The adoption of a vertical launcher greatly simplifies replenishment since each 
missile is reloaded while still inside its shipping container. 

Perhaps the single greatest advantage of the vertical launching system over 
more conventional systems is that they save valuable space--a trainable launcher 
needs clear space not only for its own rotation but also for its blast in different 
directions, which is why such launchers are so widely spaced aboard modem 
warships. In a vertical system, by way of contrast, the missile blast is concentrated 
in the immediate area of the launcher (and, to a greater extent, can be vented down 
into the launching cell), and there is no need for allowing a clear arc of fire. 
[Ref. 1] 

However, in a vertical launch surface-to-air missile (VLSAM) system, the 
aerodynamic surfaces do not suffice to control the missile at take-off, since it takes 
off at a very low initial speed. Rather, the missile generally relies on movable 
vanes set into its rocket motor exhaust. Any early failure of the booster motor will 
drop the missile directly back on its launcher. In addition, when it enters the open 
ocean environment at low velocity, the missile is subject to potentially significant 
crosswinds. The resultant of the missile and crosswind velocities is a potential high 



2 



angle of attack flow around the missile. This high angle of attack flow may cause 
asymmetric vortices to form on the missile nose and afterbody (to be discussed 
later) which induce side forces that might give flight control problems to the 
VLSAM at launch. The launch environment may contain some degree of 
turbulence, both from the atmospheric boundary layer, and from the airflow about 
the ship superstructure. Understanding the effect of this turbulence on the 
aerodynamic characteristics of a VLSAM during launch is the goal of a continuing 
effort of research and experiments conducted at the Naval Postgraduate School 
(NPS). In an earlier work, Roane [Ref. 2] developed a system of modelling 
flowfield turbulence by a series of four grids to generate turbulence in the NPS low 
speed wind tunnel. Rabang [Ref. 3] examined the effect of this turbulence on the 
asymmetric vortex induced forces generated on a VLSAM model. The goal of this 
thesis is to experimentally validate a flowfield technique to determine the location 
and intensity of the asymmetric vortices in the wake of the VLSAM model at 
various levels of freestream turbulence. The purpose is to understand better the 
asymmetric vortex behavior as effected by turbulence, and to provide flowfield 
information for correlation with numerical predictions. This study involves 
flowfield measurements at a high angle of attack about a body-only configuration 
with and without freestream turbulence. 

B. THEORY 

1. Asymmetric Vortex Theory 

The progressive development of the wake along the blunt-nosed slender 
cylindrical body when viewed in cross-flow planes is similar to the growth with 
time of the flow behind a two dimensional cylinder at an angle of attack. As shown 
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in Figure 1 from Ref. 2, a separation "bubble" exists and prevents the formation 
of the vortices in the area immediately aft the nose. Further downstream, two 
symmetrical vortices are disposed from the lee side. Theses vortices are fed by 
vortex sheets containing boundary layer fluid which has separated from the body. 
Further along the body these vortices alternately detach and move downstream one 
by one at an angle to the freestream. Other vortices form on the lee side of the body 
at an increasing distance and behave in a similar manner. This process continues 
along the body length. The result is a side force on the body, relatively small in 
magnitude, that appears as a consequence of the steady asymmetric vortices. 
[Ref. 4 & 5] 




Figure 1. Vortex Flow About an Blunt-Nosed Cylinder. [Ref. 2] 

On a slender ogive nose tip, the nose-induced separation phenomenon on a 
blunt-nosed cylinder does not occur, and consequently, nose-generated asymmetric 



4 



vortices generate much a greater proportion of the overall side force. The net 
effect of a slender body with vortex formation occurs along the entire length of the 
body, as shown in Figure 2. [Ref. 5] 




Figure 2. Vortex Flow About a Slender-Nose Cylinder [Ref. 2] 



In recent years, there has been considerable interest in the development of 
reliable techniques for the prediction of aerodynamic characteristics of missiles 
and of aircraft at high angles of attack. These research efforts have revealed that 
the flow on the leeside of these bodies is characterized by a vortex system which 
mainly depends on the angle of attack, nose shape, overall fineness ratio, crossflow 
Mach number and Reynolds number. Other secondary factors may include roll 
angle, freestream turbulence, surface roughness, acoustic environment and 
vibration. [Ref. 6, 7]. In later sections we will discuss the effects of the factors 
mentioned above. 
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a.. Vortex Regions 



As a slender body is being pitched through the angle of attack range 
0 < a < 90°, it experiences four distinct flow patterns (see Figure 3) that reflect 

the diminishing influence of the axial flow component as described by Ericsson and 
Reding [Ref. 8]. The flow regimes are characterized as follows: 

Regime I (0° < a < a sv ): At low angles of attack the axial flow 

dominates and the flow is attached. 

Regime II (Os V < a < a av ): At intermediate angles of attack, the 

crossflow sweeps the boundary layer to the leeward side where it is separated and 
rolls up into a symmetric vortex pair. The vortices are steady with time. 

Regime III (Oav < a < otu V ): At high angles of attack crossflow effects 

start to dominate and the vortices become asymmetric, thereby producing side 
forces at zero sideslip. These asymmetrical vortices gives rise to significant side 
forces and yawing moments. Nelson observed a random switching from a nearly 
symmetrical attached pair of vortices to a separated multi-vortex system occurring 
at 35° angle of attack. [Ref. 7]. The position on the slender body of asymmetrical 
vortex shedding moves forward with increased angle of attack and with Reynolds 
number and moves rearward with increasing Mach number [Ref. 6]. At the higher 
end of this angle of attack range, there are some random flow switching and flow 
instabilities. 

Regime IV ((Xuv < a < 90°): At a very high angle of attack, the 
crossflow dominates completely and the boundary layer is shed in the form of a von 
Karman vortex sheet or a random wake depending upon Reynolds number, Mach 
number and geometric details. 
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Of particular interest is Regime III where the side forces, yawing and 
rolling moments can be affected due to the asymmetrical nature of the separated 
flow field. Since these changing patterns are significant, the side forces and 
moments produced can be larger than the control moments produced by the 
deflection of the conventional control devices. And the effectiveness of aft control 
surface or fins can be greatly reduced by the vortex wake produced by forward 
separation. 

2. Two Dimensional Crossflow 

Airflow over the slender body can be divided into normal and axial 
components. The axial flow component follows along the slender body length and 
the crossflow is essentially a two dimensional flow normal to a cylinder. 

At angles of attack above 30°, in Regime III, the effective Reynolds 
number on a cylinder essentially equals the crossflow Reynolds number [Ref. 9]. 
Thus, the sectional characteristics for the slender body should be similar to those of 
a 2D cylinder. The crossflow Reynolds number is the primary factor which 
influences the separation point of the boundary layer. Viscous flow, surface 
roughness, and turbulence are other factors which influence boundary layer 
separation. Following Achenbach's terminology, a cylinder in incompressible 
crossflow experience four distinct flow regimes, each with a different type of flow 
separation, as shown in Figure 4. [Ref. 9] 

At subcritical Reynolds number, the boundary layer is laminar, and flow 
separation occurs near the lateral meridian where (p, is defined as the angle from 

the direction of the crossflow, varies 80° to 90° . 
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Figure 4. Two Dimensional Crossflow about a Cylinder [Ref. 9] 

In the critical Reynolds number range, a laminar separation bubble 
develops near the lateral meridian at (p= 90°, followed by transition in the detached 

laminar shear layer and turbulent reattachment. The reattached turbulent 
boundary layer is more energetic than the laminar boundary layer and separation is 
delayed to tp= 140°, resulting in a dramatic reduction of the drag. Thus, a drag- 

bucket is produced in the critical Reynolds number region. 

As the Reynolds number is increased, transition moves forward of the 
lateral meridian and the laminar bubble is lost. Thus, separation moves forward of 
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cp= 100° as the turbulent boundary layer thickness grows with increasing Reynolds 
number. A drag increase accompanies the wake growth. 

For asymmetric behavior in this Reynolds number regime, both 
supercritical and critical separations allow large suction peaks to be realized on one 
side of the body, whereas on the other side, subcritical separation cuts off the peak 
suction pressures, producing a large pressure differential across the body. The 
pressure drops sharply as the separation moves rearward in the critical Reynolds 
number range, resulting in a pronounced maximizing of the 2D lift to drag ratio. 
[Ref. 9] 

Figure 5 presents a logical progression of asymmetric vortex separation 
with increasing Reynolds number that explains how both a maximum and a 
minimum IC y l ma x/C n can occur in the critical Reynolds number range. At 
subcritical effective Reynolds numbers, asymmetric subcritical separation occurs 
near the 80° meridian to produce a moderate normalized side force. 

As the critical Reynolds number range is entered, critical/subcritical 
separation can occur. This provides the maximum differential in the separation 
location on opposite sides of the body and the maximum suction pressure 
differential in the vicinity of the lateral meridian where it is the most effective in 
producing a side force. The peak is sharp because a relatively small increase in 
Reynolds number results in asymmetric critical separation where nearly equal 
suction pressures occur at the lateral meridians and the separation asymmetry 
affects the pressure only at (p= 140°, where they are relatively ineffective in 

producing the side force. Thus, both the normalized side force and the side force 
itself will be very small. [Ref. 9] 
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Finally, as effective Reynolds number increases through the supercritical 
regime on to the transcritical region, the flow separation asymmetry moves toward 
the lateral meridian where it is once again efficient in generating a side force. 
[Ref. 9] 

The Reynolds number dictates the greatest influence on the normal force 
and drag characteristics, especially within the critical range where the maximum 
side force and the maximum vortex asymmetry occur. [Ref. 9] 




REYHOIDS HUMBER 



Figure 5. Side Force to Normal Force Ratio [Ref. 9] 
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3. Three Dimensional Vortices 

The nose geometry of a missile is an important factor in vortex generation 
and disposition. [Refs. 7, 10, 11, & 12] Missile nose geometry falls into two 
shapes: cones and ogives, both pointed and blunt. Noses are dimensionalized for 
comparison by their length to base diameter ratio or fineness ratio (l n /d). There 
are two kinds of asymmetric vortex shedding on slender bodies. [Ref. 13] On 
pointed bodies, the vortex asymmetry, and thus the side force, usually begins at the 
nose, and vortices are shed at a relatively rapid rate to give alternating side-force 
cells on long-slender bodies. A slight amount of nose bluntness can cause the 
second type of vortex shedding to occur. Here asymmetric vortices develop first at 
the rear of the body. and the alternating vortex shedding does not occur as 
rapidly. Thus the side force cells are much larger and can cover the entire 
cylindrical aft body. When initially generated by a slender-forebody, they are not 
affected by the aft-body vortices, though the afterbody forces will contribute to the 
overall magnitude of induced forces from the nose-body combination. 
[Refs. 7 & 9] 

For pointed conical and ogive noses, observations indicate that (Xav occurs 
as a function of the semi-vertex angle 0 a. At all Mach numbers, asymmetric vortex 

shedding starts when the angle of attack exceeds the total induced angle at the apex 
(Xav = 20 a . For a conical nose: 



0ac =6c 

For a tangent ogive nose: 



( 1 ) 
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n 



0 A = Tan" 1 ] — 

I ' 0 - 25 



( 2 ) 



For slender noses Equation 2 can be approximated by: 




(3) 



Nose fineness ratio has been shown to have an effect on the asymmetric 
vortex induced side force [Refs. 7, 10,11,12 & 14] 

The magnitude of the maximum side force generally increases with an 
increase in the nose fineness ratio and decreases with an increasing Mach number to 
nearly zero at supersonic speeds. [Refs. 15 & 16]. 

As the nose fineness ratio increases, the nose apex angle decreases and the 
angle of attack for the onset of asymmetric vortices will decrease. Thus, a missile 
with high nose fineness ratio may become more susceptible to induced side forces at 
a lower angle of attack than a missile having a lower nose fineness ratio. Jorgensen 
demonstrated that decreasing the nose fineness ratio is more beneficial in reducing 
side forces than blunting the nose [Ref. 7]. 

In addition to the nose fineness ratio, the nose-generated vortical flow is 
extremely sensitive to the nose roll angle. This phenomenon results from surface 
imperfections and deviations in the nose axisymmetric geometry. Figure 6 shows 
the work by Rabang investigating the resulting side force coefficient for nose roll 
angles varied in 45° increments. [Ref. 3] 

Should the afterbody roll angle change with nose roll angle remaining 
constant, the afterbody vortices will change the side force magnitude and direction 
as noted by Kruse [Ref. 17]. Since the nose and body of a missile do not rotate 
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independently of each other, the vortex system generated by the nose should still 
dominate afterbody vortices on an axisymmetric missile configuration. 
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Figure 6. Side Force Variations with Nose Roll Angle: Runs WOA102 

to WOA802 [Ref. 3] 

The behavior of the asymmetric vortices is well documented for a vast 
number of models and shapes, yet their cause is still not well understood. The 
general consensus among the high angle of attack researchers is that vortex 
asymmetry is induced by the location of boundary layer separation 
[Refs. 13 & 18]. But Keener and Chapman consider the principal effect to be a 
hydrodynamic (inviscid) instability in the initially symmetric vortex formation and 
the interaction of the vortices (which increases in strength with incidence) with the 
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surrounding potential flow field. In addition to the hydrodynamic instability, the 
vortex asymmetry is also strongly affected by boundary-layer (viscous) 
asymmetries resulting from transition and separation differences on each side of 
the body. [Ref. 19] 

4. Turbulence 

a. Turbulence Intensity and Length Scale 

Flow in which there are random, small-scale velocity fluctuations 
about the mean freestream flow velocity is said to be turbulent. In calculating the 
effect of turbulence on a body in the flowfield, a comparison must be made between 
the scale of the body and that of the turbulence. In addition, consideration must be 
given to the energy contained in the turbulence flowfield. 

Turbulence intensity is the measure of the relative magnitude of 
velocity fluctuations in the flowfield. Mathematically, for a horizontal flowfield or 
crosswind, the turbulence intensity T u is the ratio of the root-mean-square (rms) 
velocity fluctuation u', to the mean velocity component in the flowfield, Uoo 



T u = u’/Uoo (4) 

Since the turbulence intensity measures the magnitude of velocity 
fluctuation, a higher intensity indicates higher kinetic energy and correlates to 
more turbulent flow. 

The length scale of the turbulence, L e , is a measure of the dimension 
of the velocity fluctuation. The effect of the length scale on a body is a function of 
the relative dimensions of the turbulent eddies and the body itself. The length scale, 
despite the inference of a purely spatial measure, also represents the temporal 
character of the turbulent fluctuations. An increase in the spatial length of the 
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turbulence corresponds to an increase in the time the body is exposed to the 
fluctuation. 

Turbulence length scale to body size ratio may dictate the manner in 
which turbulence affects the VLSAM flowfield. The turbulence length scale can be 
compared to missile length, L e :Li, or missile diameter, L e : L^. When the 
turbulence length scale is many times larger than the missile length, L e » Li, the 
effect on the missile is somewhat like a uniform steady-state flowfield independent 
from the time, especially at low turbulence intensities. The effects of the flowfield 
on vortex development are largely dictated by those conditions and factors 
discussed earlier for a two dimensional cylinder. 

When the turbulence length scale is comparable to the body length, 
L e = Li, the flowfield is distinctly non-uniform and may cause unwanted rolling, 
pitching and yawing motion. [Ref. 20] In both cases, the missile body motion 
results from the actual flowfield conditions about the missile body. 

When the turbulence length scale is smaller than the body, L e « Li, 
the result has a distinctly different effect. Most importantly, when it is comparable 
or smaller than the missile diameter, L e < Lj, the small scale turbulence present in 
the flowfield may change the boundary layer and flow separation present on the 
missile surface, particularly when the turbulence length scale is the same size or 
smaller than the boundary layer thickness. The results of Rabang show an increase 
in turbulence intensity with length scale on the order of the nose-generated vortex 
scale tends to delay the onset of induced side force, to steady the asymmetric vortex 
formation, and to leave the maximum value C y more or less unaltered. [Ref. 3]. 
The results of Rabang show an increase in turbulence intensity with a length scale 
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on the order of the boundary-layer scale tends to reduce the magnitude of the 
induced side forces. [Ref. 3] 

Both large and small scale turbulence length scales actually coexist in 
the flowfield since the turbulence length is the mean of the disturbances in the 
flowfield. The combination of different length scales may add complexity to the 
qualitative analysis of the flowfield. As the disturbances in the flowfield experience 
strain, the disturbances will break up into finer disturbances of smaller scale and 
decreasing energy. The process repeats itself or cascades until the small 
disturbances eventually disappear owing to viscosity. When the disturbances 
decrease in size and energy, their individual intensities will decrease at a faster rate. 
Thus, the turbulence length scale average will be biased towards the larger 
turbulence scales in the flowfield and gives an accurate quantitative representation 
of all the length scales within the flowfield, since the energy and intensity of the 
larger scale turbulence is dominant [Ref. 21.]. 




Figure 7. Vortex Shedding on a Missile at High Angle of Attack 

[Ref. 18.]. 
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5. Effects of Body, Wing, Strakes and Tail 

The vortex structure of the complete missile can be constructed by 
superimposing those of the body, wing strakes, and tail. This results in a vortex 
system as shown in Figure 7. [Ref. 18] Body vortices interact with wings and tails, 
and wing or canard vortices interact with tails. The interaction can be very strong 
if a vortex passes very close to a tail surface. The addition of wings causes the nose 
vortices to move closer to the body. As the vortices move closer to the body, the 
result is comparable to increasing the effective angle of attack causing unsteady 
asymmetric vortices. The net effect of the wing-body combination seems to be a 
reduction in the effective angle of attack for the onset of asymmetric vortices and 
side forces. [Ref. 22] 

For the body-only case, from the previous research by Rabang, an 
increase in turbulent intensity and scale tend to reduce the induced side forces and 
to shift the range of action to slightly higher inside angles. [Ref. 3] 

The addition of strakes will generate additional vortices, improving the 
interactive between body vortices and airflow over the wings. The placement of 
long strakes along the missile body and nose was found by some researchers to 
reduce the forces and moments generated by asymmetric vortices. [Ref. 18]. 
Rabang demonstrated the addition of wings and strakes typical of a VLSAM tended 
to preserve the induced side force for all levels of turbulence intensities and length 
scales. [Ref. 3] 

The adding of a tail has, in general, relatively little influence on the 
maximum side forces. With the application of wings of high aspect ratio, the tail 
loses its efficiency at a small angle of attack, but gains in efficiency with the 
growing angle of attack. [Ref. 18] 
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The use of small nose bluntness, nose booms and boundary layer trips has 
proven effective in reducing asymmetric vortex-induced loads. [Ref. 23] 

C. EFFECT ON VERTICALLY-LAUNCHED SURFACE TO AIR 
MISSILE 

1. Marine Environment 

Conditions in the marine atmospheric environment may have a significant 
effect on the VLSAM. The atmospheric boundary layer (ABL) is the result of the 
interaction of the atmospheric flow over the underlying land or sea surface. This 
layer is characterized by a turbulent transfer of momentum, heat and mass (water 
vapor) and their associated gradients. 

The surface layer, approximately the lowest 10% of the A.B.L., is 
characterized by the mechanically produced turbulence from the surface roughness 
and friction with nearly constant vertical fluxes of momentum, heat and mass. The 
height of this layer may vary over a wide range but is typically of the order of 
50 meters. In the lower most part of the surface layer, commonly referred to as 
the roughness layer, the surface has the greatest effect on generating turbulence and 
influencing fluid motion. Mean flow in the roughness layer is nonhomogeneous 
and three-dimensional in nature. In the marine roughness layer, the traveling 
surface waves affect the flow which result in various properties of the sea surface. 
However, the observation of these wave effects has been hindered by the motion of 
the instrument platform induced by the wind and waves as well as by the presence 
of the platforms (towers, buoys, ships) in the flow field. Furthermore, the 
majority of the flow in the surface layer itself can be considered horizontally 
homogeneous. [Ref. 20] 
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A measure of the general roughness of the surface is given by Zq, the 
roughness length. Its value is determined as a function of the mean wind velocity at 
various elevations above the surface of interest. Surface layer turbulence length 
scale and intensity can be empirically determined by combining the roughness 
length with the elevation and windspeed. [Ref. 24] 

For a typical open ocean roughness length in the range 10' 3 < Z Q < 10‘ 2 
meters, turbulence intensities are on the order of 13 to 17 percent at a 10 meter 
elevation and mean wind speed of approximately 25 m/sec. [Ref. 24:p. 1 1] 

Apparently, the turbulence intensity in the marine surface layer can be 
significant. However, the effect of such turbulence fluctuations is highly dependent 
upon the length scales present. For the conditions stated in the preceding 
paragraph, the longitudinal turbulence length scales are in the range 80 < L e < 90 
meters. For a typical missile with a 1.1 foot diameter, the turbulence length scale to 
missile body scale ratio is about L e : = 280:1, and would have little effect on its 

boundary layer development. Yet, the cascade effect presented above allows for a 
decrease in the length scale until viscous forces dissipate the energy. Therefore, it 
is possible that the initial large scale turbulence and crosswind interaction with the 
ship superstructure would decrease, or "cascade," to scales where they could affect 
the development of the missile boundary layer. 

2. Launch and Crosswind Velocities. 

A typical VLSAM launching with a lOg acceleration reaches an altitude of 
approximately 56 feet traveling at a vertical velocity of V v = 164 ft/sec. The 
VLSAM is still well within the surface layer environment (50 meters) and is 
subject to crosswind and turbulence effects. As an example, if a ship travels at 20 
knots and a mean wind speed of 20 m/sec are combined with the VLSAM launch 
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velocity, the resultant vector is at 191 ft/sec at an angle of 31° from the missile 
flight direction. The result represents the VLSAM flying at an effective angle of 
attack of 31° at 191 ft/sec. This places the VLSAM in Regime III, the asymmetric 
vortex region, as it travels through the surface layer. [Ref. 2] 

Gregoriou also states that an effective angle of attack up to a = 50° may be 

reached by a VLSAM while it pitches over towards a target. [Ref. 25] These 
examples illustrate the distinct possibility of asymmetric vortex induced side forces 
on a VLSAM during launch phase. 

3. Other Launch Considerations 

There are many factors which affect the aerodynamics of VLSAM during 
the launch phase of its flight. They can be divided into two categories. One 
category of factors are inherent to the design of the missile, such as the plume (or 
jet) effect of the missile's engine, blast effects of the vented exhaust gases and 
activation of the flight control system. The other category of factors includes 
shipboard roll, pitch, and yaw, and ship airwake turbulence. [Ref. 2] 

Plume and blast effects can affect the VLSAM while the propulsion 
system gases or exhaust gases impinge on the aerodynamic surfaces of the missile. 
In addition, the exhaust gases can directly impact the flowfield into which the 
missile is launched, especially if the gases are vented upward into the vicinity of the 
accelerating missile. 

The VLSAM itself may impose some form of flight control, aerodynamic 
or thrust vector, to maintain flight attitude during the launch phase. [Ref. 25] 
Should the missile change flight attitude, the flowfield around it will be altered. 

The ship undergoes constant roll, pitch and yaw motion and transmits this 
motion to the launch platform of VLSAM— the canister, since the canister used for 
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storage is also the launcher. This increases the complexity to the initial velocity 
vector as the VLSAM leaves the guide. 

In addition to the turbulence generated in the marine atmospheric 
boundary layer, turbulence may be produced by the ship's hull and superstructure. 
This flowfield or ship airwake may increase turbulence intensities while decreasing 
turbulence length scales. 

The above considerations, while not forming a complete description, 
emphasize the numerous factors which may influence the flight behavior of the 
VLSAM during the launch phase. Except for the limited discussion in the 
preceding paragraphs, these additional effects are beyond the scope of the present 
investigation. 
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II. EXPERIMENT APPARATUS 



A. EXPERIMENTAL HARDWARE 

The experimental hardware for this experiment can be divided into 7 major 
subdivisions. There are the wind tunnel, the VLSAM model with associated 
support structures, the 3-D traverser, the scanivalve mechanism, the turbulence 
generating grids, the 5-hole pressure probe and the data acquisition system. Each 
of these will be individually discussed in the following sections. 

1. Wind Tunnel 

The horizontal-flow, low-speed wind tunnel that was used during this 
research is located at the Naval Postgraduate School and is shown in Figure 8. 
[Ref. 26:pp. 3-7] 

It is a single return tunnel measuring 64 feet in length and between 21.5 
and 25.5 feet in width. The power section of the tunnel comprises a hundred 
horsepower electric motor coupled to a three-blade variable pitch fan by a four- 
speed Dodge truck transmission. Immediately following the propeller blades are a 
set of stator blades to help straighten the flow. Turning vanes are installed at all 
four comers, with two small-mesh wire screens installed upstream of the settling 
chambers, to smooth the airflow and reduce turbulence. A heavy wire screen is 
also installed downstream of the test section to prevent damage to the fan blades or 
turning vanes in the event of model failure during operation. The test section of the 
tunnel measures 45 inches by 32 inches, with comer fillets which house the test 
section lighting; the cross sectional area is 9.88 square feet. The contraction ratio 
of settling chamber area to test section area is approximately ten. The walls of the 
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test section diverge slightly to allow for boundary layer growth without a reduction 




Figure 8. Naval Postgraduate School Low Speed Wind Tunnel 

[Ref. 26:p. 3-7] 

As the test section is designed to operate at atmospheric pressure, a 
circumferential breather slot is installed downstream of the test section to replenish 
air lost through leaks in the tunnel walls. The tunnel was designed to provide 
velocity up to 290 feet per second in the test section. [Ref. 26] 

A reflection plane is installed in the floor of the test section, which decreases the 
vertical dimension to 28 inches. In the center of the reflection plane is a flush- 
mounted turnable for adjustments in pitch angle. The angle of the turnable is 
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remotely controlled with an electric motor installed beneath the tunnel. 
[Ref. 26:pp. 3-3, 3-4] See Figure 9. 




Figure 9. Turnable Assembly and VLSAM in Wind Tunnel 
Flow measurement components utilized with the wind tunnel are a dial 
thermometer, manometer, and pitot static tube. Temperature in the wind tunnel is 
measured by a dial thermometer extending into the setting chamber. Dynamic 
pressure in the wind tunnel is measured by the calibrated static pressure difference 
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between the test section and the settling chamber using a water filled manometer. 
Static pressure in the settling chamber is measured by four static pressure taps, one 
on each wall , upstream of the test section to preclude interference from the model. 
The pressure taps at each section are connected via a common manifold prior to 
feeding into the manometer. The manometer measures pressure differences in 
centimeters of water, which is the test section dynamic pressure, and is converted to 
the actual wind tunnel velocity by Equation (5). 



A pressure transducer circuit connected to the manometer displays test 
section dynamic pressure on a digital readout. 

2. VLSAM and Support Equipment 

The VLSAM model was constructed of 6061 and 2024 aluminum alloy by 
Naval Postgraduate School personnel. It is intended to model a cruciform tail- 
control missile with very low aspect ratio long dorsal fins. It is designed to permit 
force and moment measurements while operating in a subsonic wind tunnel. 



! New calibration factor. 




(5) 



where 



• U m = measure velocity (ft/sec) 

• 2.0475 = conversion factor. 

• Pcm H 2 O = manometer reading in cm of H 2 0 

• K = Calibration factor 1 , 0.889 for no grid, 1.5084 for 



Grid #1. 

• p = air density(lb/ft 3 ) 
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The body section is a hollow cylinder with locating pin attachment points 
for the balance gage, sleeve, wings and tails. The machined sleeve provides a close 
tolerance fit between the balance gage and the interior of the VLSAM model . The 
model is seated to the balance gage by locating pins. Body roll angle can be varied 
in 45° increments. Nose roll angle may be varied in 45° increments independent of 
the body. The nose is a pointed ogive shape. Four very low aspect ratio wings with 
strakes comprise the cruciform wing section and four tail control fins are mounted 
equilaterally in fixed position along the model axis. All parts are rigidly connected 
to the model body by countersunk machine screws. The surface of the VLSAM 
model is buffed to a polished finish and is free from protuberances. The 
dimensions of the model are summarized below as shown in Figure 10. 

The VLSAM model support is rigidly fixed in the wind tunnel test section 
by the reflection plane turntable at the base and an aluminum reinforced clear 
plexiglass section at the top. The plexiglass sheet has slots cut to the width of 5/4 
inches and the length of 7, 8, and 10 inches.These slots correspond to the positions 
of length/diameter ratios of 3, 6, and 9 of the missile model, i.e., 5.25, 10.5 and 
15.75 inches from the nose, respectively. 

A planar survey grid, the Y-Z plane, is perpendicular to the freestream 
velocity and down stream 10.5 inches from the nose of the missile. The magnitude 
of the Z dimension is divided evenly with respect to the center line of the missile 
body (see Figure 11). The Y-Z dimensions for the preliminary run and the actual 
runs are 6.5 x 8 and 3 x 5.5 inches respectively. The step distance s are 0.5 inch and 
0.25 inch for the preliminary run and for the actual runs respectively. Because of 
the probe measuring volume, the grid mesh could not be finer than 0.25 inch. 
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Total length 
Base diameter 
Length/diameter ratio 
Ogive nose length 
Ogive/diameter ratio 
Wing span/root chord 
T ail span/root chord 
Center of pressure 



22.85 in. 

1.75 in. 

13.06 
4.0 in. 

2.29 

3.13 in./ 1 3.55 in. 

5.50 in. /1. 70 in. 

13.5 inches all of nose tip (approx.) 



NAVAL POSTGRADUATE SCHOOL 
SURFACE-TO-AIR 
MISSILE MODEL 





Figure 10. Drawing of VLSAM Model and the Specification [Ref. 3] 
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(b) Front View 

Figure 11. The Planar Survey Grid. 



29 




Figure 12. VLSAM Model with Support Equipment and Pressure 
Probe in the Test Section (without Wing) 

The sting model support changes model pitch angle via a chain gear drive 
powered by an electric motor in the horizontal direction with the pivot point 
coinciding with the approximate center of the VLSAM model (see Figure 12). 

3. Velmex 8300 3-D Traverser 

The Velmex 8300 is a 3-D traversing system using three microcomputer- 
controlled stepping motors (one for each axis of movement). The system is 
composed of the motor controller assembly and the traversing assembly as shown 
in Figure 13(a),(b). 

A 5-hole pressure probe was attached to the 8300 control drive. It can 
accurately and effectively move the pressure probe through the test section. The 
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motor controller assembly is the interface unit between the operator and the 
motors. The controller is capable of interpreting motor movement commands 
from a host computer, programmable control, or terminal. Data communication 
between a terminal and the controller assembly is made through a full-duplex RS- 
232C port located on the front of the controller assembly. Manual motor 
movements are input via switches also located on the front panel. 

The resident BASIC interpreter software program contains the necessary 
motor movement subroutines and is responsible for monitoring motor status. Fault 




Figure 13(a). Velmex 8300 3-D Traverser - Motor Controller 

Assembly 
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Figure 13(b). Velmex 8300 3-D Traverser -- Traversing Assembly 
lights on the front panel indicate parity or data bus errors and a ready light 
indicates the previous motor movement is completed and the controller is ready to 
accept the next input. Additionally, the controller has 12k of custom firmware for 
interactive motor control. The combined 20k of ROM is contained on a 6502 based 
microcomputer along with 4k of RAM. Operator-selectable motor variables, 
accessible through software commands, include motor velocity, motor 
acceleration, increment distance, and incremental unit (motor steps or inches) 
[Ref. 27]. 

The traverser assembly is a stainless steel and aluminum assembly 
consisting of three separate motor / jackscrew assemblies. Each motor/jackscrew 
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moves the traverser along one of the axes depending on the connections made to the 
controller. Each motor is a 200 step per revolution, 10 amp stepper motor with a 
maximum velocity of 3000 steps/sec. Minimum motor movement is one motor step 
(1/200 of a revolution) which equates to approximately 0.000125 inch. [Ref. 27] 

An extension aluminum bar was attached to the base of the traverser and 
mounted to existing hardware located on top of the tunnel. The assembly was 
mounted in a position such that tunnel-induced vibrations on the traverser were 
minimized. Additionally, a thin vertical extension bar was inserted behind the 
probe tubing to reduce the pressure probe vibrations due to the high dynamic 
loading in the high velocity airstream. 

Motor movement were controlled via the controller assembly from an 
IBM PC/AT computer located adjacent to the tunnel. The program for traverser 
movement will be discussed in detail in the software section. 

4. Scanivalve 

One 48-port scanivalve was used to measure each of the 5-hole pressure 
probe pressures. (See Figure 14.) The scanivalve was used because it provides an 
economical means of sensing the pressures of a multiple number of sources. It also 
simplifies calibration procedures in that only one transducer needs to be calibrated. 

The scanivalve mechanism puts out a 7-bit binary coded decimal (BCD) 
signal that corresponds to the port (1-48) currently connected to the scanivalve 
transducer. This arrangement allows remote, electronic monitoring of the 
configuration of the port assembly of the scanivalve. 

The scanivalve consists of two separate units. One unit contains the 
transducer, motor drive, and port assembly. The second unit houses the solenoid 
controller. Control of the scanivalve is accomplished through the solenoid 
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controller. The solenoid controller allows two commands, STEP and HOME. 
STEP "steps" the scanivalve one port location. HOME command sends the 
scanivalve to port number 48. (See Figure 15.) 



Figure 14. 48-Port Scanivalve 

Computer control is accomplished through the solenoid controller by 
means of the Hewlett-Packard (HP) Relay Actuator. It will accept a 15 voltage 
pulse to either STOP or HOME. 
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Scaniualue Control 





Figure 15. Scanivalve Control 

The signal path for the scanivalve, containing information about the 
5-hole probe port pressure, passes through the signal conditioner/amplifier to the 
Relay Multiplexer; then the values measured by the HP Digital Multiplexer are sent 
to the computer. The signal is conditioned by a low pass filter with a cut off of 
10Hz. This removes most of the noise, including power line interference, prior to 
being measured of the HP Digital Multiplexer. The transducer signal is fed through 
a 1000 gain low noise amplifier to improve resolution for small magnitude signals. 
The HP Digital Multimeter converts the analog voltage signal into a digital signal 
for use by the data acquisition program named "PPROBE." An IBM-AT 
microcomputer executes the the PPROBE software and stores the collected data in 
the consecutive files. 

5. Turbulence-Generating Grids 

The turbulence grids were designed for use in the low speed wind tunnel 
to generate turbulence in varying intensities and length scales [Ref. 2]. (See 
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Figure 16). Each grid was mounted in a wood frame, which was placed 73 inches 
ahead of the pivot axis of the model support system (see Figures 17). 




Figure 16. Turbulence Grids 2, 3, and 4, Clockwise from the Left 

Grid 2, 3, and 4 
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Figure 17. Planview of VLSAM Model and Pressure Probe in the Test 

Section (not to scale) 



TABLE 1. GRID SPECIFICATIONS 1 



Grid 


Mesh Width (in.) 


Bar Diameter (in.) 


Mesh/Diameter 


Material 


One 


5.00 


1.00 


5 


Wood 


Two 


3.75 


0.75 


5 


Wood 


Three 


2.50 


0.50 


5 


Wood 


Four 


1.00 


0.0625 


16 


Wire 



The specifications for the four grids are listed in Table 1 and Figure 1 8. 
Three of the four grids are square-mesh square-bar and fabricated from wood. 



![Ref. 2:p. 28] 
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They employ bi-plane construction which generates nearly isotropic homogeneous 
turbulence. The fourth grid is also square-mesh, but was fabricated from round 
wire. Turbulence intensities and length scales for each of the grids were 
determined by Roane [Ref. 2] and are shown in Figures 1 9 and 20. 




_L 

T 



Bar Diameter (d) 



Me3h 

Width 



Figure 18. 



Square-Mesh Turbulence-Generating Grid 
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2 The dashed line indicates the model pivot axis. [Ref. 2:p. 45] 
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Figure 20. Grid Turbulence Length Scales 3 
The dynamic pressures shown in Table 2 are the actual dynamic pressures 
in the test section from Roane [Ref. 2, p. 44-49]. The conclusion from Rabang 
[Ref. 3] was that the effects of the grid generated turbulence with regards to 
changing length scales at constant intensity or changing intensities with constant 
length scale could not be investigated with the present grid turbulence parameters. 



3 The dashed line indicates the model pivot axis [Ref. 2:p. 48] 
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TABLE 2. GRID TURBULENCE PARAMETERS: AT THE MODEL 

PIVOT AXIS 4 



Grid 


Intensity 

(percent) 


Length Scale (in.) 


Turbulence/ 
Model Dia. 


Dynamic 

Pressure 

(lb./ft 2 ) 


One 


3.31 


1.84 


1.05 


15.35 


Two 


2.78 


1.56 


0.89 


14.88 


Three 


1.88 


1.08 


0.62 


16.38 


Four 


0.47 


0.27 


0.15 


15.61 


None 


0.23 


- 


- 


15.85 



6. 5-Hole Pressure Probe 



Since the local flow angles are high just outside of the vortex core, and the 
3-sensor hot wire probe is limited to flow angularities of no more than 45°, from 
the results of experiments by Naik it was shown that a pressure probe is useful to 
measure the flowfield where the flow angularity is high (typically greater than 35°) 
[Ref. 28]. The data from the pressure probe can be reduced to obtain the 
distributions of local total pressure, static pressure and velocity vectors 
simultaneously. The wide measurement range is the chief advantage the pressure 
probe has over other measuring devices. A 5-hole pressure probe is used in this 
experiment. The disadvantage is the limitation of time-averaged measurements. 

The 5-hole pressure probe is made of corrosion resistant non-magnetic 
stainless steel silver bragged. It is 0.125 inch in diameter and 24 inches in total 
length with 22 inches of reinforcement tubing. There are five take-off tubes with a 
reinforcement block on the top. At the measuring tip is a five hole prism shaped 



4 [Ref. 2:pp. 44-49] 
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measuring section. (See Figure 21) The pressure probe gives the data of pitch 
angle, yaw angle and velocity. More details can be found in Ref. 29. The pressure 
data was reduced to obtain isobars of total pressure coefficient and static pressure 
coefficient, and to map the crossflow velocity vectors. 




Figure 21. 5-Hole Probe and its Measuring Tip [Ref. 29] 
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The speed of reading depends on the length and diameter of the pressure 
passage inside the probe, the size of the pressure tubes to the manometer, and the 
displacement volume of the manometer. The time constant increases rapidly for 
smaller diameter tubes. The diameter of tube used in this experiment is 1/4 inch 
O.D. and the length of the tubes are three feet, so the time delay is about 0.15-0.26 
second. 

7. HP Data Acquisition System 

The HP data acquisition system consists of a combination of hardware and 
software that allows the IBM PC-AT computer to act as a fully automated 
instrumentation system. [Ref. 29] The individual HP instruments used in this thesis 
include the Relay Multiplexer, Digital Multimeter and Relay Actuator. (See 
Figure 22.) Each of the instruments can be operated manually at the computer 
screen by means of the mouse control. 

In this mode the instrument’s operating controls and functions, as well as 
digital displays, are relayed to the computer screen. Instead of user interaction 
with the instrument's controls and indicators, each unit is manipulated by the 
computer mouse controller. 

In addition to the manual mode, the HP data acquisition system can be 
operated in the program (Basic) mode. This method is employed in this thesis. A 
program was written in the advance Basic language (BASICA) that called up each 
instruments and its function as necessary. All data acquisition related to the 
VLSAM wind tunnel experiment was accomplished by the use of the BASICA 
program. 

The Relay Multiplexer provides one common output channel which in 
turn can be read by the Digital Multimeter. 
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Figure 22. HP Data Acquisition System 
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Digital Multimeter (DMM) Soft Front Panel. [Ref. 30] 
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The Digital Multimeter (DMM) can measure + or - DC voltages, AC 
voltages and Ohms. In this thesis only the DC voltage measurement function is 
used. An auto range is used in which the DMM selects the optimal range for the 
signal that is being measured. The DMM automatically converts input analog 
voltage signals into a digital (or binary) form which can be read by the computer. 
The DMM has a continuous data sampling rate of 2.5 or 12.5 readings per second 
and the higher sampling rate, 12.5 readings per second, with a accuracy of ±0.05% 
of the input voltage, is used throughout the experiments. (See Figure 23.) 

The Relay Actuator is used solely in controlling the scanivalve to STEP or 

HOME. 

The whole setting combined with the 3D-traverser, scanivalve, 5-hole 
pressure probe, HP data acquisition system and IBM PC AT is shown in Figure 24. 

B. EXPERIMENT SOFTWARE 

Several programs were used for data acquisition, data reduction and plotting in 
this experiment. The relation between those programs is shown in Figure 25. Each 
of them is discussed as follows. 

1. PPROBE Program 

The BASICA application program that runs the VLSAM experiment is 
composed of three parts. 

a. ST ATE. FILE 

The first is called a STATE.FILE. The statefile is a program 
automatically compiled by the HP Soft Front Panel software. It tells the computer 
what configuration each unit was left in when the instruments were last used (i.e., 
Relay Multiplexer was set for channel one as the input and the output device was 
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enabled etc.) The statefile in use in the VLSAM laboratory is named 
PGMSHEL.HPC and is located in the Lung sub-directory on the computer's fixed 
C-disk. [Ref. 30] 

Data Requisition S u stem 




Figure 24. Whole Hardware Setting 
b. PGMSHEL 

The second part of the BASICA application program is called 
PGMSHEL. PGMSHEL consists of BASICA program lines that perform 
initialization chores to allow communication between the HP instruments and the 
IBM computer. In essence, the PGMSHEL allows the computer to know every 
function available at each of the data acquisition instruments. When one of these 
functions is called up in the BASICA application program the computer already 
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"knows" that function exists and where to find it. PGMSHEL, like the 
STATE.FILE, is, created by the HP system soft Front Panel Software.[Ref. 30]. 




Pro g ram Procedures 



Figure 25. Program Procedures 
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c. Traverser 



The third part of the BASICA application program is the actual 
application code. PGMSHEL occupies lines 1 through 999. The application code 
starts at line 1000 and begins to run only after all system initialization is completed. 
The application code consists of a combination of HP instruments statements with 
BASICA keywords, the scanivalve control program and the traverser program. 

The traverser program was originally written by Kindelspire 
[Ref. 31] in Advanced Basic (BASICA) language to serve as the interface between 
the operator and the motor controller assembly. It allows the operator to control 
precisely the pressure probe movement via the 3-D traverser controller unit. It 
consists of manual input (meaning one motor movement for each operator input) 
and computer-controlled movement. It was minimally modified by the author in 
manual input and the author developed a new computer-controlled movement. The 
algorithm flow chart is shown in Figure 26 for this experiment. The PPROBE 
program is shown in Appendix A. 

The manual movements were used to pre-position the pressure probe 
prior to the data collection run through the whole measured plane. 

The operator first uses manual input to move the total pressure hole 
(Pi) of the pressure probe as near as possible in the horizontal direction to the 
center axis of the VLSAM model, then he moves the probe vertical down to the 
position where he desires to measure. This point is called the original point, which 
is used for the reference at each run. After pre-positioning the pressure probe to 
the original point, the operator selects the computer-controlled movement option. 
At this option, the operator is prompted to input the dimension measurement in 
(Y,Z) format and the step distance. All units are in inches. The program will 
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automatically calculate and display the number of points in the Y-direction and 
Z-direction and the total number of points to be measured. If the operation is 
satisfied, then the operator is prompted to input the name of the data file to be 
stored; otherwise, the program goes back to prompt the operator to re-input the 
dimension measurement. After the operator inputs the name of the data file, the 
program displays the name of the data file which automatically increases by the 
point number in the Y-direction with a extension, DAT. One column of data is 
written to a single file. 5 The program prompts the operator to have a final check on 
the dimension measurement, the number of points and the name of the data file 
before the data can be acquired. Before the data acquisition, the operator must 
rotate the pressure probe until the pressure for the ports P2 and P3 are equal 
(nulling); then he inputs the yaw angle indicated by the traverser scale wheel. The 
first step was to move the scanivalve from port 1 to port 4 6 via the Relay Actuator 
with a one second delay time which permits the pressures to equalize before the 
DDM samples the output voltage from the scanivalve transducer via the Relay 
Multiplexer. 

After ten samples are taken, the Relay Actuator steps the scanivalve to 
the next port and the process is repeated until all the five channel pressures are 



5 If there are total 16x 10 = 160 points to be measured and the input the name of the 
data file is TEST, then the date file will automatically increase from TEST 01. DAT 
to TEST 16. DAT through the whole measurement. 



6 Scanivalve port 4 is consistent with Pi of the pressure probe and port 5 is 
consistent with P2 and so on. 
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PPROBE Pro g ram Floii> Chart 




Figure 26 PPROBE Program Flow Chart 

measured (i.e., from port 4 to port 8). Then the Relay Actuator homes the 
scani valve to port 48. The program will show all the readings and the average 
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values for each channel. The program prompts the operator to check whether the 
data is valid or not. If not, it will remeasure; otherwise, the program will move the 
pressure probe upward one step automatically. It then repeats the data acquisition 
procedure. After the program finishes the one column data acquisition, it stores all 
the averaged data of each sample point in a data file, then moves the traverser to the 
next column position. The program repeats the data acquisition process until all the 
columns are finished. 

2. CALP Program 

Before and after each time the CALP program is executed, the transducer 
voltage output is calibrated by applying a known pressure source. From the output 
voltage and pressure data, the slope equation is calculated by curve fitting.The slope 
equation is used by the PPROBE program to convert the transducer voltage to 
physical pressure. The CALP program is shown in Appendix B. 

3. CONVERT Program 

The program is written by the author in Fortran. It opens a consecutive 
file created by program PPROBE and reads the contained data in each row and then 
converts the data into the Y-Z coordinates, velocity, yaw angle, pitch angle, total 
local pressure, the total local pressure coefficient, local static pressure and the local 
static pressure coefficient which is stored into a file named "RESULT.DAT" for 
later use by the REVX program. The CONVERT program is shown in Appendix C 
and the pressure probe calibration chart is shown in Appendix D. The data 
reduction can be divided into 5 parts for interpretation and the flow chart is shown 
in Figure 27. 
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Conuert Flou> Chart (T) 
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P i : initial etompheric pressure (in Hg) 

p 

f : final atompheric pressure (in Hg) 

T i : initial temperature (°F) 

Tj : final temperature (°F) 

K : wind tunnel calibration factor 



^atm : atompheric pressure (psf) 

Tave : average temperature (°R) 
P : air density (lb/ft 3 ) 



P=Delta P(Y) *2.0475 + P. 

atm 










Delta P(Y): convert voltage data to 
physical pressure data 








Calculate the pitch 
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FPitch(P): pitch angle function | 







FYSLOP : velocity function | 










-40°< 0C<- 1 0° 


-iO # iaiio* io*<<x<40* 














Calculate local 
velocity 




YAW : yaw angle 

FYAW : yaw angle function 



Figure 27(a) CONVERT Program Flow Chart 
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Conuert Flom Chart (2) 




Q : Local Dynamic Pressure 

Q t : Free Stream Dynamic 
Pressure 



FYa^Jtev^Ang^ 



P TC =FPT(cc) 



Calculate the local total 
pressure coefficient 



Calculate the local 
static pressure 
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FPT : Total Pressure Function 



P TC : Pitch Angle Pressure 
Coefficient 

P T : Local Total Pressure 



Write No, X, Y, Yel, Yaw, 
pitch, P T , CP t , P s , CP S 
to RESULT.DAT 



C End ^ 



Pj : Local Total Pressure 

CP t : Local Total Pressure 
Coefficient 

Ps : Local Static Pressure 
CP S : Local Static Pressure 
Coefficient 



Figure 27(b) CONVERT Program Flow Chart 
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a. The Coordinate 

The Y-Z coordinate is the same as in the original files. 

b. The Pitch Angle 

The pitch angle coefficient 0c is determined by calculating Equation 



, P4-P5 

c "Pi-P 2 



( 6 ) 



where 0c = pitch angle coefficient 

• Pi = indicated total pressure 

• P2 = indicated static pressure 

• P4 = pitch angle pressure 

• P5 = pitch angle pressure 

After the pitch angle coefficient is obtained, the corresponding pitch 



angle must be found. It is calculated by using Equation 7 which was developed by 



the author by using commercial curve-fitting software to fit the provided 



calibration chart. 



0 = 3.759 + 53.75680c - 1.3O850 C 2 - 1.65830 c 3 - O.8O610 C 4 - 16.51 150 c 5 (7) 

where 

• 0 = pitch angle (degree) 

• 0 C = pitch angle coefficient 

c. Velocity: 

At any particular pitch angle, 0, the velocity pressure coefficient can 
be read from the velocity pressure coefficient curve on the calibration chart. The 
velocity pressure coefficient curve can be represented by Equations 8, 9, and 10 at 
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different pitch angle regions. They were also developed by the author by using 
commercial curve-fitting software. 

At0 <-10° 



V c = 0.981-0.01020 - 3.0*10-402 - 2.5*10-603 (8) 

At-lO°<0<lO° 

V c = 0.98-0.0060 +2.0*10-402 (9) 

At0 >10° 



V c = 0.981-0.00350- 1.143*10-402 - 5.833*1O 6 03 (10) 

where 

• 0 = pitch angle (degree) 

• V c = velocity coefficient 

And the velocity pressure coefficient is represented by Equation 1 1 . 



V c 



Pt-Ps 



£ 




P1-P2 P1-P2 P1-P2 



where Pt = total pressure (psf) 

• P s = static pressure (psf) 

• p = gas density ((lb/ft 3 ) 

• Pi = indicated pressure (psf) 

• P 2 = indicated static pressure (psf) 
. V = velocity (ft/sec) 

• q = local dynamic pressure 



Recall that 



( 11 ) 



Patm — pRT 
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atm 



p “ R*T 

where 

• Patm = ambient pressure (psf) 

♦ R = gas constant =1718 (ft-lb/slug*R) 

♦ T = temperature (°R) 

• p = density (lb/ft 3 ) 

Combine Equations 11 and 12. The velocity is calculated by Equation 13 



( 12 ) 



v=y 



2V c (P,-P 2 ) 



(13) 



and adding the wind tunnel calibration factor, K, and the conversion factor gives 
the local velocity in the wind tunnel: 



u m =Y 



2*2.0475*V c *(Pi-P2) 

p*K 



(14) 



d. Pressure Coefficient 

At a particular pitch angle, 0, the total pressure coefficient, P tc , can 
also be read from the total pressure coefficient calibrate curve. At different pitch 
angle regions, it can be represented by Equation 15. 



Ptc - 0. 1 



0< -30° 

-30°<e<-20° P Xc = 0.02 + 1*10' 3 0 



At S 



-20° < 0 < 30 e 



(15) 



P tc =° 



30° < 0 < 40° ? Xc = 0.03 - 1*1 O’ 3 0 



where 
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• Ptc = the total pressure coefficient of pressure probe calibration 

• 0= pitch angle 
and 



PrPt 



PtL- 




where 

• Pi = indicated total pressure 

• PtL = Local total pressure 

• P s l = Local static pressure 

PtL' PsL = jpU m 2 



( 16 ) 



(17) 



where U m is from Equation 14. From Equations 16 and 17, both P^and P s l are 
calculated. 



Since the pressure measurements are affected by slight shifts in wind 
tunnel velocity and temperature, pressure coefficients are used when the local and 
total tunnel pressures are compared and they are non-dimensionalized by the tunnel 
dynamic pressure. These coefficient are defined below : 

Cp.= £ ^r 1 (is) 

Cps = ■ Ps Q Ps (19) 

where 

• Q = Freestream dynamic pressure 

• Pt = Freestream total pressure 

• P s = Freestream static pressure 

• Cpt = Total pressure coefficient. 

• Cps = Static pressure coefficient. 
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e. Yaw Angle 

Yaw angle is simply read from the traverse scale wheel when the 
pressures of the probe tip were nulled. This value is input to the program. 

/. Temperature 

Temperature is taken from the average of initial and final 
temperature. Runs are limited in time to allow for a maximum temperature 
excursion of 15°F. 



Tav ~" 



max 



min 



g. Pressure 

Pressure is read from the barometer and was recorded before and 
after each run. 

4 . REVX Program 

This program is used to reverse the Y-coordinate of RESULT.DAT and 
to store the data in PLOT.DAT for the XPLANE program. The REVX program is 
shown in Appendix E. 

5. XPLANE Program 

The purpose of this program is to read the data of PLOT.DAT and 
convert those data into the crossflow direction and magnitude for the ARROW 
program. The XPLANE program is shown in Appendix F. 

6. ARROW Program 

This program plots the crossplane velocity vector plots. The ARROW 
program is shown in Appendix G. 



58 



C. EXPERIMENTAL CONDITIONS 



To correlate the data with those of Rabang, all the experimental conditions 
were kept as much like Rabang's as experimental conditions would allow. 

Reynolds number, test section dynamic pressure, test section velocity and nose 
geometry and roll angle are kept constant as possible. 

Dynamic pressures were set to match those treated by Rabang, giving a 
subcritical Reynolds number of R e = l.lxlO 5 . 

The angle of attack was fixed at 50° to get comparative data with that of Rabang 
for which maximum side force occurred at the no wing and no grid conditions. (See 
Figure 28.) 




Figure 28. ROB802: Side Force Coefficient [Ref. 3] 
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For all cases of Rabang, the rose roll angle was held fixed. Afterbody roll 

angles of the reference system as shown in Figure 29 [Ref. 3]. 

• Body A: wings and tail at zero roll angle in a "+" configuration. 

• Body B: no wings or tails, body roll angle set at <J>r = 45°. 

• Body C: wings and tails at (()r = 45° roll angle in a "x" configuration. 






Figure 29. Body Configurations and Reference System [Ref. 3] 

The longitudinal position for measurements was at a length/body diameter 
ratio of 6. This position equates to 10.5 inches from the nose of the missile model. 

The goal of this experiment was to validate the pressure measurement system 
for flowfield variables with elevated levels of turbulence; to determine the location 
and intensity of the asymmetric vortices in the wake of the VLSAM model at a 
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raised level of freestream turbulence; and to display the asymmetric vortices by 
velocity mapping and pressure contours. The purpose is to correlate the results 
with the force measurements of Rabang to provide a greater understanding of the 
vortex flowfield. 

D. EXPERIMENT PROCEDURE 

1. Dynamic Pressure Calibration 

For the purpose of continuing research, all of the turbulence grids were 
calibrated in the tunnel. Normally, the pressure difference between a set of static 
ports in the settling chamber and a set of static ports in the test section entrance is 
used to relate to the actual dynamic pressure in the test section. Because of the 
disturbance on the second set of static ports caused by the frames which support the 
grids, each grid must be tested for its own tunnel calibration constant. Readings 
from the tunnel calibration manometer and from a pitot-static tube inserted in the 
center of the test section were recorded over a speed range, and a calibration curve 
was generated for each grid. The calibration curves will be discussed in the Results 
section. 

2. Scanivalve Transducer Calibration 

The transducer voltage is first adjusted about zero in in milli voltage 
range, then the "CALP" program is run. 

The scanivalve is calibrated by applying a known pressure source to the 
scanivalve transducer which is connected to the HP Digital Multimeter instrument. 
The known pressure source was provided by the calibration manometer. (See 
Figure 30.) The examination value of the calibration pressure is displayed on the 
Digital Panel Meter (DPM) (see Figure 31) and the readings are input to the 
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computer manually at the same time. Both transducers are referenced to an 
ambient pressure. A HP Digital Multimeter converts the analog voltage signal 
from the scanivalve transducer into a digital signal and stores the collected data into 
a calibration file. This procedure is performed before and after each run to 
measure any shifts between the initial and final readings. 

These data were used to generate the calibration curve equation which 
relates the transducer voltage output to the wind tunnel dynamic pressure. 

The calibration curve equation is used in the data reduction computation 
to transform voltage data into dynamic pressure data. 
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Figure 30. Calibration Manometer 
3. Pressure Probe Calibration 

The pressure probe was calibrated by the manufacturer prior to use in this 
experiment. The calibration chart is shown in Appendix D. The curves of velocity 
pressure coefficient, pitch angle, pressure coefficient and total pressure coefficient 
vs pitch angle have been transformed to polynomial equations by the author. This 
curve-fitting has been discussed in CONVERT program section. The 5-hole 
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pressure probe is rotated until opposing peripheral ports measure equal pressures. 
The corresponding angle rotation determines the yaw flow angle. For this reason, 
it is referred to as a nulling probe. Although the nulling procedure allows for 
analysis of large yaw angle flow, the procedure is mechanically complex and time 
consuming. 




Figure 31. Digital Panel Meter (DPM) 



4. Preliminary Run and Data Collection 

The purpose of the preliminary run was to determine the location of the 
asymmetric vortices. The dimension of 6.5 by 8 inches and the step of 0.5 inch and 
the no grid condition were used in the preliminary run. Prior to each run, the 
initial settling chamber temperature and barometric pressure were recorded and 
the data acquisition system and traverser were fully energized and allowed to settle 
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out for about half an hour. During the settling period, the data acquisition system 
was cycled repeatedly by the program to bring the relay actuator, relay multiplex, 
digital multimeter, scanivalve amplifier and traverser to the operating temperature 
condition and to check the voltage output until all readings were consistent. Next, 
the wind tunnel was started and the flow was set at a different velocity for a small 
test measuring probe and the recorded data were checked for consistency. If they 
were not constant, the data acquisition system was re-calibrated again. The 
condition of each run is shown in Appendix H . The time of the preliminary run 
was about 6 hours. In order to keep the air density of test section as uniform as 
possible, the experimental procedure was interrupted when the settling chamber 
temperature rose 15°F above the operating temperature. The experiment was 
continued when the settling temperature had cooled down within 5°F of the 
operating temperature by the exchange of air. At the end of each run, the final 
settling chamber temperature and the barometric pressure were recorded again and 
the average value of them was used in the CONVERT program. 

5. Actual Run 

From the preliminary result, the dimension of 3 by 5.5 inches and a step 
of 0.25 inch were used through the experiment. This dimension covers the main 
portion of asymmetric vortices. Because the experiment requires considerable 
time for the nulling of the probe tip pressures before each measurement can be 
taken, and because the settling chamber temperature tended to rise quickly due the 
friction of the air flow through the turbulence generating grid, the time duration 
was about 10 ~1 1 hours for each run. 



65 



E. EXPERIMENTAL CORRECTIONS 



Blockage factor corrections were made by Equations 20 and 21 from Ref. 26. 



u = U„(1+E) 


(20) 


q = q u (l+2e) 


(21) 



where 

• U: test section velocity (lb/ft 2 ) 

• U u : Measured test section velocity (lb/ft 2 ) 

• q: dynamic pressure (lb/ft 2 ) 

• q u : Measured dynamic pressure (lb/ft 2 ) 

• e: blockage correction 

The total blockage correction can be expressed as Equation 22. 

1 model frontal area .... 

£ = T : (22) 

4 test section area v 

The blockage equation for each body configuration as a function of the model 

angle of attack has been investigated by Rabang as shown in Figure 32. The angle 

of attack was fixed for the experiments at 50°. The value of e equals 0.0123. These 

equations are implemented into the data conversion program listed in Appendix B. 
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ANGLE OF ATTACK (DEGREES) 



Figure 32. Blockage Factors [Ref. 3] 
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III. RESULTS 



A DYNAMIC PRESSURE CALIBRATION 

Wind tunnel dynamic pressure calibration was conducted prior to the 
experiment for each grid and no grid. From these data, the wind tunnel calibration 
factors w'ere obtained. This reference dynamic pressure is used for adjusting the 
tunnel flow velocity to the expected experimental condition for the different grids. 
The calibration curves and equations for each run are shown in Table 3 in Fig. 33. 



TABLE 3. WIND FUNNEL CALIBRATION RESULT 



No Grid 


Grid #1 


Grid #2 


Grid #3 


Grid #4 


Ptrue Ph:0 


Ptrue Ptuo 


Ptrue Phuo 


Ptrue 


PpbO 


Ptrue 


Pf-bO 


0 


0 


0 


0 


0 


0 


0 


0.36 


0 


0.64 


0.56 


0.61 


1.17 


0.8 


0.75 


0.66 


0.89 


0.99 


0.19 


0.94 


2.11 


3.34 


2.63 


1.75 


1.68 


0.98 


2.14 


1.75 


1.26 


1.33 


3.83 


4.29 


3.82 


2.6 


2.55 


1.55 


3.47 


2.63 


2.75 


3.33 


5.93 


6.65 


4.77 


3.12 


3.62 


2.20 


4.97 


3.58 


4.12 


4.62 


7.94 


8.84 


7.56 


5.02 


4.55 


2.82 


6.38 


4.45 


5.52 


5.90 


9.97 


11.25 


9.80 


6.53 


5.62 


3.43 


8.13 


5.50 


6.84 


7.07 


11.40 


12.80 


11.93 


7.91 


6.74 


4.09 


9.25 


6.05 


9.16 


9.10 










8.13 


4.95 


10.08 


6.45 


11.0 


10.37 










9.19 


5.66 


















10.04 


6.17 
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Figure 33(a). Dynamic Pressure Calibration at no Grid 




Figure 33(b). Dynamic Pressure Calibration at Grid #1 
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Figure 33(c). Dynamic Pressure Calibration at Grid #2 




Figure 33(d). Dynamic Pressure Calibration for Grid #3 
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Figure 33(e). Dynamic Pressure Calibration at Grid #4 
Figures 33 show the tunnel calibration factors k are 0.8891, 1.5084, 1.6487, 
1.6545 and 1.1 167 for no grid, grid 1, grid 2, grid 3 and grid 4 respectively. Due 
to the particular portable manometer, the curves were shifted to go through the 
origin. These values were used in the CONVERT program for calculating the 
pressure and velocity in the test section. 



B. TRANSDUCER CALIBRATION 

Transducer calibration was conducted both before and after each run. The 
results of the two calibrations for each run are presented in Table 4 and Figure 34 
which shows a slight shift between initial and final runs. For a given voltage, the 
initial and final equations differed by approximately 5.16% for no grid and 0.03% 
for grid 1 for an experiment time duration of about 10 hours. An average of the 
two equations was used in the CONVERT program for all case runs to minimize the 
error. In each run, when the voltage difference between P 2 and P 3 is within 0.05, 
the data were recorded. The maximum relative error is 1.7% in pressure and 
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0.85% in velocity for Grid 1. For no grid condition, the maximum relative error 
is 3.23% in pressure and 1.5984% in velocity. 

For this experiment, since the air flow Mach number was 0.1, the 
incompressible flow assumption was a reasonable approximation. In addition, it is 
assumed that the temperature and velocity gradients are negligible in a small region 
around the probe. 



TABLE 4. TRANSDUCER CALIBRATION RESULT 





TEST RUN 


G#0 


G#1 




pressure 


voltage 


pressure 


voltace 


pressure 


voltage 




-8.91 


-2.06663 


-8.26 


-1.9361 


-8.22 


-1.9152 




-4.37 


-1.03319 


-3.44 


-0.8304 


-3.92 


-0.9304 


Before 


0.29 


0.08565 


0.55 


0.1278 


1.22 


0.3017 




4.87 


1.17548 


3.24 


0.7727 


5.64 


1.3487 




8.80 


2.1055 


6.18 


1.4682 


10.39 


2.4762 








9.92 


2.3514 


13.52 


3.2192 




-8.81 


-2.00795 


-9.26 


-2.1754 


-7.77 


-1.8257 




-4.46 


-1.01094 


-4.37 


-1.0658 


-3.53 


-0.8461 


After 


0.49 


0.172168 


-0.19 


-0.0720 


4.89 


1.166 




4.74 


1.17583 


6.41 


1.5074 


9.16 


2.1786 




9.52 


2.3122 


9.87 


2.3242 


12.32 


2.9292 
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Voltage (Test run) 



yl= ■ 02472 + 42352x 

0 Pressure(BT) 

• Pressure(AT) 

y2 = - 0.0898 + 4.23 19x 




y = 0.0214 + 4.2321 x 
B 0 ’i+y2)/2 



Figure 34(a). Transducer Calibration at no Grid Condition for the 

Preliminary Run 
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Data from "Transducer Calib. Data" 



yl= -0.0115 + 4.2277X 

D Pressure(BO) 

♦ Pressure(AO) 

y2 = 0.0543 + 4.2365x 



-3 -2 - 1 0 1 2 3 

Voltage (G^-0) 



y = - 0.1685 + 4.2335X 

0 (yi+yW 



0 2 4 6 8 10 12 

Voltage (G#=0) 

Figure 34(b). Transducer Calibration at no Grid Condition 
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Data from "Transducer Calib. Data" 

yl= - 0.0635 + 4.2244 x 

° Pressure(Bl) 

♦ Pressure(AJ) 

y2 = -0.0186 + 4J161x 

-2-1 0 1 2 3 4 

Voltage (G#=l) 








y = - 0.041 + 4.2203x 
D (yi+y2)!2 



Figure 34(c). Transducer Calibration at Grid #1 Condition 



C. PRELIMINARY RUN 

The objective of the preliminary run was to determine the location of 
asymmetric vortices. This result was used to determine the size of the measured 
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plane and the step distance for the following runs, so that more accurate velocity 
and pressure data could be obtained. The preliminary run conducted 6d 
downstream of the nose of the missile in the plane normal to the freestream at 
ambient turbulence conditions, consisted of a 6.5- by 8-inch measured phase, a step 
distance of 0.5 inch, and the no grid condition. Figure 35 shows the velocity 
vectors superimposed on the measured grid which clearly shows the location of the 
vortex. The nose induced asymmetric vortices are marked by the swirling vectors. 
The velocities give a representation of an indirect flow visualization of the vortex. 
The crossflow is expected to reach zero at the freestream region. This 
inconsistency in arrow magnitude and direction are seen at the outer boundaries 
because of a possible misalignment of the traverser. Figures 36 and 37 show the 
2D contour and 3D surface plots of local total pressure coefficient and local static 
pressure coefficient with the same data. The value between two contour lines is 
0.15. From Figure 36, it is clear that the location of asymmetric vortices are 2.5 to 
5.5 inches on the Z axis and 0 to 2 inches on the Y axis. Hence a small planar survey 
grid, 3 x 5.5 inches and with a step of 0.25 inch, was determined for further 
experiments. 

Figures 36 and 37 show the total pressure coefficient to be uniform at the 
freestream and a high negative pressure gradient to exist in the core of the vortex. 
The static pressure coefficient has a positive value at the center of the measure 
planar and has a moderate negative pressure gradient in the vortex core. More 
accurate investigations in the vortex region are discussed in the following section. 
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Figure 35. Preliminary Crossflow Velocity Vector Superimposed on 

the Grid. 
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Y AXIS (IN) 



TEST RUN, TOTAL PRESSURE COEFF. 




Figure 36(a). Total Pressure Coefficient Contour 7 



7 Dashed line indicates the negative value. Note: data as plotted times 100. The 
distance between two isobars is 15. 
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TEST RUN, STATIC PRESSURE COEFF. 




Figure 36(b). Local Static Pressure Coefficient Contour 8 



8 Dashed line indicates the negative value. Note: data as plotted times 100. The 
distance between two isobars is 15. 
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TOTAL PRESSURE COEFF. 



TEST RUN, NO GRID 




Figure 37 (a). 3D Surface Plot. Local Total Pressure Coefficient 
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TEST RUN, NO GRID 




Figure 37 (b). 3D Surface Flot. Local Static Pressure Coefficient 

D. BODY ONLY WITHOUT TURBULENCE 

The results from Rabang show a maximum side force coefficient to be about - 
4.2 at 50° AOA for body-only and the no grid condition. [Ref. 3] (See Figure 28.) 

As marked by the swirling vectors in Figure 38, two vortices form two circles 
and rotate in opposite directions. Represented by the velocity magnitude, the 
vortex strength reaches its maximum in the outer core of the two vortices. The 

i 
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maximum side force is due to the relative strengths of the two vortices and the 
displaced location of the vortices from a symmetric configuration. 

The total pressure coefficient is expected to be constant in the outer 
incompressible and inviscid portion of the vortex. The viscous forces found within 
the vortex core decrease the flow’s fluid-mechanical energy; therefore the total 
pressure coefficient becomes negative. The vortex contour with a C pt of -1.35 is 
approximately 0.57d above the missile surface. The other vortex with a C pl of 
-1.43 is approximately 0.26d above the missile surface. It is clear to see from 
Figure 39(a) that the extent of each vortex is about 0.86d wide and the distance 
between vortex centers is about 0.86d. In a negative sense, closer isobars have a 
higher gradient which forms two inverted spike shaped vortices as shown in Figure 
40(a). It is clear to see the shapes of vortex cores in this figure. 

In the freestream condition, the static pressure coefficient is expected to be 
zero. [Ref. 32] While the static pressure coefficient progresses inward from the 
freestream conditions to the vortex core, the angular velocity increases, static 
pressure drops, and the static pressure coefficient itself decreases from zero to a 
negative value. The vortex contour with a C ps of -1.8 vortex is approximately 
0.37d above the missile surface. The other vortex contour with a C ps of -2.22 is 
approximately 0.29d above the missile surface. Figure 39(b) shows 1.1 4d for the 
width of each vortex. The static pressure coefficient has a lower gradient inward to 
the vortex core. Figure 40(b) exhibits a fatter inverted spike shape, indicating a 
diffused and lesser strength. The vortex closer to the missile is stronger, as can be 
seen by the higher local gradients and the deeper peak. 
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Figure 38. Crossflow Velocity Vector Superimposed on the Grid for 

no Grid Condition 
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Y AXIS (IN) 



NO GRID, TOTAL PRE55URE COEFF. 




Figure 39(a). Total Pressure Coefficient Contour at no Grid 

Condition. 9 



9 Dashed line indicates the negative value. Note: data as plotted times 100 The 
distance between two isobars is 15. 
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NO GRID, STATIC PRESSURE COEFF 




Figure 39(b). Local Static Pressure Coefficient Contour at no Grid 

Condition. 10 



10 Dashed line indicates the negative value. Note: data as plotted times 100. The 
distance between two isobars is 15. 
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TOTAL PRESSURE COEFF. 



NO GRID 




Figure 40(a). 3D Surface Plot. Local Total Pressure Coefficient at no 

Grid Condition. 
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NO GRID 




Figure 40(b). 3D Surface Plot. Local Static Pressure Coefficient at no 

Grid Condition. 



E. BODY-ONLY WITH TURBULENCE 

Rabang demonstrated that for the body-only case, an increase in turbulence 
intensity tends to reduce dramatically the maximum induced side force. [Ref. 3] 
(see Figure 41). 
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Figure 41 R1B801. Side Force Coefficient. [Ref. 3] 

Figure 42 shows the crossflow velocity vectors and a slight shift in direction. 
Represented by the swirling vectors, the vortex has less strength than in the no 
turbulent condition. The asymmetry can still be seen in Figure 42. 

The total pressure coefficient has a fatter spike shape and more of a plateau at 
the vortex core than in the no turbulent condition as shown in Figure 44(a). One 
vortex has a C pt contour of -0.97 a distance of 0.43d above the missile surface and 
the other one has a C pt contour of -1.16 a distance of 0.4d above the missile surface. 
The distance between the two vortices' cores is 0.7 Id, which is 0.1 5d closer than the 
one in the no turbulent condition. (See Figure 43(a).) Each vortex has a larger 
extent than the one with no turbulence at the same isobars. 

The static pressure coefficient, similar to the total pressure coefficient, has a 
more diffused shape than in the no turbulent condition. (See Figure 4(b).) The C ps 
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contour of -0.45 for one vortex is more or less 0.36d above the missile surface and 
the C ps contour of -0.72 for the other vortex is approximately 0.5d above the 
missile surface. The distance between these two vortices’ cores is closer than that in 
no turbulent condition; each vortex has a width of about 1.1 4d. (See Figure 43(b).) 

Comparing the extents of the vortices with and without flowfield turbulence, a 
more tightly wrapped vortex (smaller diameter) has a greater strength. The result 
of the turbulence is to diffuse, but not to eliminate, the asymmetric vortices. With 
the additional turbulence, the vortex centers are closer to the missile surface and 
more symmetric. These facts agree with Rabang’s results. 
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Figure 42. Crossflow Velocity Vector Superimposed on the Grid for 

Grid 1 Condition. 
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Y AXIS (IN) 



GRID=1, TOTAL PRESSURE COEFF. 




Figure 43(a). Total Pressure Coefficient Contour at Grid 1 

Condition. 11 



n The dashed line indicates the negative value. Note: data as plotted times 100. The 
distance between two isobars is 15. 
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Y AXIS (IN) 



GRID= 1 , STATIC PRESSURE COEFF. 




Figure 43(b). Local Static Pressure Coefficient Contour at Grid 1 

Condition 12 



12 Dashed line indicates the negative value. Note: data as plotted times 100. The 
distance between two isobars is 15. 
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TOTAL PRESSURE COEFF. 



GRID=1 




Figure 44(a). 3D Surface Plot. Local Total Pressure Coefficient at 

Grid 1 Condition. 
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STATIC PRESSURE COEFF. 



GRID=1 




Figure 44(b). 3D Surface Plot. Local Static Pressure Coefficient at 

Grid 1 Condition. 
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IV. CONCLUSION AND RECOMMENDATIONS 



The flowfield downstream of a vertically-launched surface-to-air missile 
model at an angle of attack of 50° and a Reynolds number of 1.1 x 10 5 was 
investigated in a wind tunnel of the Naval Postgraduate School. The body-only 
configuration was tested. Two flowfield conditions were treated: the nominal 
ambient wind tunnel condition, and a condition with grid- generated turbulence of 
3.8% turbulence intensity and a dissipation length scale of 1.7 inches. The 
following conclusions were reached: 

1) The relative strengths of the asymmetric vortices can be noted by the sharp 
spike shape in the ambient condition; this condition becomes diffused and becomes 
fatter in the turbulent condition; 

2) The right side vortex has greater strength than the left side one as seen by 
the changes in the total pressure coefficient and static pressure coefficient contours 
with and without a turbulent condition; 

3) An increase in turbulence intensity tends to reduce the strength of the 
asymmetric nose-generated vortices; 

4) An increase in turbulence intensity also pushes two asymmetric vortices 
closer together; 

5) Crossflow velocities were examined and were found to indicate the 
behavior denoted by the pressure contours. 

To understand better the behavior of asymmetric vortices induced by the 
missile nose, a few recommendations are suggested as follows: 

1 ) Examine the vortices at positions of length/diameter ratios of 3 and 9; 
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2) Investigate the interaction effects of wings and consider the effect of roll 
angle on asymmetric vortices; 

3) Investigate the asymmetric vortex behavior affected by various angles of 
attack and turbulent conditions with varying intensities and length scales. 
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APPENDIX A. PPROBE PROGRAM 



1 DEF SEG:CLEAR ,&HFEOO:GOTO 4 'BEGIN PCIB PROGRAM SHELL 

2 GOTO 1000 ’ USER PROGRAM 

3 GOTO 900 ' ERROR HANDLING 

4 I=&HFE00 ’ COPYRIGHT HEWLETT-PACKARD 1984,1985 

5 PCIB.DIR$=ENVIRONS("PCIB") 

6 IS=PCIB.DIRS4'VCIBILC.BLD" 

7 BLOAD IS,I 

8 CALL I(PCIB.DIRS,I%,J%):PCIB.SEG=I% 

9 IF J%=0 THEN GOTO 1 3 

10 PRINT "UNABLE TO LOAD,"; 

1 1 PRINT " (ERROR #";J%;")" 

12 END 

13 ’ 

14 DEF SEG=PCIB.SEG:O.S=5:C.S=10:I.V=15 

15 I.C=20:L.P=25:LD.FILE=30 

1 6 GET.MEM=35:L.S=40:PANELS=45:DEF.ERR=50 

17 PCIB.ERR$=STRINGS(64,32) : PCI B.NAMES=STRINGS( 16,32) 

18 CALL DEF.ERR(PCIB.ERR,PCIB.ERRS,PCIB.NAMES,PCIB.GLBERR) : PCIB.BASERR=255 

19 ON ERROR GOTO 3 

20 J=-l 

21 IS=PCIB.DIRS+"\PCIB.SYN" 

22 CALL O.S(IS) 

23 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

24 1=0 

25 CALL I.V(I,READ.REGISTER,READ.SELFID,DEF1NE,INITIALIZE.SYSTEM) 

26 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

27 CALL I.V(I£NABLE.SYSTEM,DISABLE.SYSTEM,INITIALIZE,POWER.ON) 

28 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

29 CALL I.V(I,MEASURE,OUTPUT,START,HALT) 

30 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

31.CALLI.V(I r ENABLE.INT.TRIGGER,DISABLE.INT.TRIGGER£NABLE.OUTPUT,DISABLE.OUTT) 

32 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

33 CALL I.V(I,CHECK.DONE,GET.STATUS,SET.FUNCTION,SET.RANGE) 

34 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

35 CALL I.V(I,SET.MODE,WRITE.CAL,READ.CAL,STORE.CAL) 

36 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

37 CALL I.V(I,DELAY,SAVE.SYSTEMJJ) 

38 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

39 1=1 

40 CALL I.V(I,SET.GATETIME,SET.SAMPLES,SET.SLOPE,SET.SOURCE) 

41 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

42 CALL I.C(ITREQUENCY,AUTO.FREQJ>ERIOD,AUTO.PER) 

43 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

44 CALL I.CaJNTERVAL,RATIO,TOTALIZE3100MILLI) 

45 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

46 CALL I.C(I,R1,R10,R100,R1KILO) 

47 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

48 CALL I.C(I,R10MEGA,R100MEGA,CHAN.A,CHAN.B) 

49 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

50 CALL I.C(I, POSITIVE, NEGATIVE.COMN.SEPARATE) 
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51 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

52 1=2 

53 1=3 

54 CALL I.V(I,ZERO.OHMS,SET.SPEED,J,J) 

55 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

56 CALL I.C(I,DCVOLTS .ACVOLTS ,OHMS .R200MILLI) 

57 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

58 CALL I.C(I,R2,R20,R200,R2KILO) 

59 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

60 CALL I.Ca£20KILO,R200KILO,R2MEGA,R20MEGA) 

61 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

62 CALL I.C(I,AUTOM,R2.5,R12.5,J) 

63 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
641=4 

65 CALL I.V(I,SET.COMPLEMENT,SET.DRIVER, OUTPUT.NO. WAIT, ENABLE.HANDSHAKE) 

66 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

67 CALL I.V(IT>ISABLE.HANDSHAKE,SET.THRESHOLD,SET.START.BIT,SET.NUM.BITS) 

68 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

69 CALL I.VCLSET.LOGIC.SENSEJJJ) 

70 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

71 CALL I. C(I .POSITIVE .NEG ATI VE.TWOS.UNSIGNED) 

72 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

73 CALL I.C(I, OC.TTL.R0, Rl) 

74 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

75 CALL I.C(I,R2,R3,R4,R5) 

76 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

77 CALL I.C(I,R6,R7,R8,R9) 

78 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

79 CALL I.C(I,R10,R1 1.R12.R13) 

80 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

81 CALL I.C(I,R14,R15,R16,J) 

82 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

83 1=6 

84 CALL I.V(I,SET.FREQUENCY,SET.AMPLITUDE,SET.OFFSET,SET.SYMMETRY) 

85 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

86 CALL I. V(I,SET.BURST.COUNT J J J) 

87 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

88 CALL I.C(I,SINE,SQUARE,TRIANGLE,CONTINUOUS) 

89 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

90 CALL I.C(I,GATED,BURST4 J) 

91 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 
921=7 

93 CALL I.V(I,AUTOSCALE,CALIBRATE,SET.SENSmVITY,SET.VERT.OFFSET) 

94 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

95 CALL I.V(I,SET.COUPLING,SET.POLARrrY,SET.SWEEPSPEED,SET.DELAY) 

96 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

97 CALL I. V(I, SET .TRIG .SOURCE, SET .TRIG.SLOPE.SET.TRIG.LEVEL.SET.TRIG.MODE) 

98 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

99 CALL I.V(I,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO,GET.TIMEBASE.INFO) 

100 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

101 CALL I. V(I,GET. TRIG. INFO .CALC.WFVOLT .CALC. WFTIME.CALC.WF.STATS) 

102 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

103 CALL I.V(I,CALC.RISETIME,CALC.FALLTIME,CALC.PERIOD,CALC.FREQUENCY) 

104 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

105 CALLI.V(I,CALC.PLUSWIDTH,CALC.MINUSWIDTH,CALC.OVERSHOOT,CALC.PRESHOOT) 
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106 IF PCIBERRoO THEN ERROR PCIB.BASERR 

107 CALL I.V(I,CALCEK.TO.PK,SET.TIMEOlJT,SCOPE.START,MEASURE.SINGLE.WF) 

108 IF PCIBERRoO THEN ERROR PCIB.BASERR 

109 CALL I.Va.MEASURE.TWO.WFJJJ) 

1 10 IF PCIBERRoO THEN ERROR PCIB.BASERR 

111 CALL i.c(ieionano,rioonanoeimicroeiomicro) 

112 IF PCIBERRoO THEN ERROR PCIB.BASERR 

113 CALL I.C(IE100MICROE1MILLLR10MILLIE100MILLI) 

114 IF PCIBERRoO THEN ERROR PCIB.BASERR 

115 CALL I.Ca^l^lOE20NANOE200NANO) 

1 16 IF PCIBERRoO THEN ERROR PCIB.BASERR 

117 CALL I.C(I,R2MICRO,R20MICRO,R200MICRO,R2MILLI) 

118 IF PCIBERRoO THEN ERROR PCIB.BASERR 

119 CALL I.CaE20MILLIE200MILLI,R2E20) 

120 IF PCIBERRoO THEN ERROR PCIB.BASERR 

121 CALL I.Ca,R50NANOE500NANOE5MICROE50MICRO) 

122 IF PCIBERRoO THEN ERROR PCIB.BASERR 

123 CALL I.C(IE 500MICRO ,R 5MILLI ,R 50MILLI E500MILLI) 

124 IF PCIBERRoO THEN ERROR PCIB.BASERR 

125 CALL I.C(I t R5,R50,CHAN.A,CHAN.B) 

126 IF PCIBERRoO THEN ERROR PCIB.BASERR 

127 CALL I.C(IEXTERNALEOSITIVE,NEGATIVE,AC) 

128 IF PCIBERRoO THEN ERROR PCIB.BASERR 

129 CALL I.C(I.DC, TRIGGERED, AUTO.TRIG.AUTO.LEVEL) 

130 IF PCIBERRoO THEN ERROR PCIB.BASERR 

131 CALL I.C(I,X1.X10,STANDARD,AVERAGE) 

132 IF PCIBERRoO THEN ERROR PCIB.BASERR 

133 1=8 

134 CALL I. V(I,OPEN .CHANNEL, CLOSE.CHANNELJJ) 

135 IF PCIBERRoO THEN ERROR PCIB.BASERR 

136 CALL C.S 

137 IF PCIBERRoO THEN ERROR PCIB.BASERR 

138 IS=PCIB .DIRS+"\PCIB .PLD” 

139 CALL L.P(IS) 

140 IF PCIBERRoO THEN ERROR PCIB.BASERR 

141 IS="DMM.01":I=3:J=0:K=0:L=1 

142 CALL DEFINE(DMM.01,IS,IJ,KE) 

143 IF PCIBERRoO THEN ERROR PCIB.BASERR 

144 IS=’’FUNC.GEN.0r:I=6:J=0:K=l:L=l 

145 CALL DEFINE (FUNC GEN.0 1 JS,U JC E) 

146 IF PCIBERRoO THEN ERROR PCIB.BASERR 

147 I$="SCOPE.01 ":I=7:J=0:K=2:L= 1 

148 CALL DEFINE(SCOPE.01,I$4d,KE) 

149 IF PCIBERRoO THEN ERROR PCIB.BASERR 

150 I$="COUNTER.01":I=1:J=0:K=3:L=1 

151 CALL DEFINE(COUNTER.01 ,ISJ J,KE) 

152 IF PCIBERRoO THEN ERROR PCIB.BASERR 

153 I$="DIG.IN.01":I=4:J=0:K=4:L=1 

154 CALL DEFINE(DIG.IN.01,I$,U4CE) 

155 IF PCIBERRoO THEN ERROR PCIB.BASERR 

156 I$="DIG.OUT.0 1 " :I=4:J= 1 :K=4:L= 1 

157 CALL DEFDME(DIG.OUT.0U$,U,KE) 

158 IF PCIBERRoO THEN ERROR PCIB.BASERR 

159 I$="RELAY.ACT.01":I=8:J=0:K=5:L=1 

160 CALL DEFINE(RELAY.ACT.0U$,UEE) 
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161 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

162 IS="RELAY.MUX.01":I=2:J=0:K=6:L=1 

163 CALL DEFINE(RELAY.MUX.01,I$,I,J,K,L) 

164 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

800 I$=ENVIRONS("PANELS")+'\P ANELS.EXE” 

801 CALL L.S(IS) 

899 GOTO 2 

900 IF ERR=PCIB.BASERR THEN GOTO 903 

901 PRINT "BASIC ERROR #" ;ERR;" OCCURRED IN LINE ";ERL 

902 STOP 

903 TMPERR=PCIB .ERR :IF TMPERR=0 THEN TMPERR=PCIB.GLBERR 

904 PRINT "PC INSTRUMENT ERROR #”;TMPERR;" DETECTED AT LINE ";ERL 

905 PRINT "ERROR: ";PCIB.ERRS 

906 IF LEFTS(PCIB.NAMES,l)oCHRS(32) THEN PRINT "INSTRUMENT: ”;PCIB.NAMES 

907 STOP 

908 COMMON PCIB.DIRS.PCIB.SEG 

909 COMMON LD.FILE,GET.MEM,PANELS,DEF£RR 

910 COMMON PCIB.BASERR,PCIB.ERR,PCIB.ERRS,PCIB.NAME$,PCIB.GLBERR 

911 COMMON READ .REGISTER .READ.SELFID, DEFINE, INITIALIZE.SYSTEM, 

ENABLE.SYSTEM,DISABLE.SYSTEM,INmALIZE,POWER.ON,MEASURE, OUTPUT, 

START,HALT£NABLE.INT.TRIGGER,DISABLE.INT.TRIGGER£NABLE.OUTPUT, 

DISABLE.OUTPUT,CHECK.DONE,GET.STATUS 

912 COMMON SET.FUNCTION,SET.RANGE,SET.MODE,WRITE.CAL£EAD.CAL, 

STORE.CAL,DELAY,SAVE.SYSTEM,SET.GATETIME,SET.SAMPLES,SET.SLOPE, 

set.sourcezero.ohms,set.speed,set.complement,set.driver, 

OLTPUT.NO. WAn\ENABLE.HANDSHAKE,DISABLE.HANDSHAKE 

913 COMMON SET .THRESHOLD, SET.START.BIT,SET.NUM.BITS,SET.LOGIC.SENSE. 

SET.FREQUENCY,SET.AMPLITUDE,SET.OFFSET,SET.SYMMETRY, 

SET.BURST .COUNT, AUTO SCALE, CALIBRATE,SET.SENSITIVITY, 
SET.VERT.OFFSET,SET.COUPLING,SET.POLARITY,SET.SWEEPSPEED 

9 14 COMMON SET.DELAY,SET.TRIG.SOURCE,SET.TRIG.SLOPE,SET.TRIG.LEVEL, 

SET.TRIG. MODE, GET.SINGLE.WF, GET .TWO.WF.GET.VERT.INFO, 

GET.TIMEBASE.INFO,GET.TRIG.INFO,CALC.WFVOLT,CALC.WFTIME, 

CALC.WF.STATS,CALC.RISETIME,CALC.FALLTIME,CALC.PERIOD 

915 COMMON CALCJREQUENCY, CALC.PLUS WIDTH, CALC.MINUSWIDTH, 

C ALC.O VERS HOOT,C ALC PRES HOOT.C ALC.PK.TO.PK.SET.TIMEOUT, 

SCOPE.ST ART, MEASURE.SINGLE.WF.MEASURE.TWO.WF, OPEN .CHANNEL, 
CLOSE.CHANNEL 

916 COMMON FREQUENCY,AUTO.FREQ,PERIOD,AUTOEER,INTERVAL, RATIO, 

TOTALIZER 100MILLI.R 1 £10£ 100.R 1KILOE IOMEGA^ 100MEG A, CHAN .A, 
CHAN.BPOSm VE.NEGATIVE.COMN, SEP AR ATE, DC VOLTS, AC VOLTS, OHMS, 
R200MILLI£2£20JI2(X)£2KILOJI20KILO£200KILO 

917 COMMON R2MEG A.R20MEG A,AUTOM,R2.5,R 1 2.5 ,POS ITIVENIEG ATIVE.TW OS , 

UNSIGNED ,OC,TTL,R0,Rl,R2,R3,R4,R5,R6,R7,R8,R9,R10,Rll,R12,R13,R 14, R 15, 
R16,SINE,SQUARE, TRIANGLE, CONTINUOUS, GATED, BURST.R10NANO, 

R 100NANOJU MICROS 10MICROE 100MICRO 

918 COMMON R1MILLI£10MILLI£100MILLI£1£10£20NANO£200NANO£2MICRO, 

R20MICRO,R200MICRO,R2MILLI,R20MILLI,R200MILLI,R2,R20,R50NANO, 
R500NANO£5MICRO£50MICRO£500MICRO£5MILLI£50MILLI£500MILLI£5, 
R50 .CHAN. A.CHAN.BJEXTERNAL POSITIVE 

919 COMMON NEGATIVE.AC, DC, TRIGG ERED,AUTO.TRIG,AUTO.LEVEL,X 1 10, 

STANDARD AVERAGE 

920 COMMON DMM.01,FUNC.GEN.01,SCOPE.01,COUNTER.01 £>IG.IN.01,DIG.OUT.01, 
RELAY. ACT .01 41ELAY.MUX.01 

999 END PCIB PROGRAM SHELL 
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1000 REM THIS STEP INITIALIZES THE HP SYSTEM 
1010 CLS 

1020 OPTION BASE 1 

1030 DIM P(5)PA(50,5),PP(50,5),XPT(50),YPT(50),X(50),Y(50),YAW(50) 

1040 REM 

1050 CALL INITIALIZE.SYSTEM(PGMSHEL.HPC) 

1060 REM 

1070 REM SET FUNCTIONON THE ’DMM' , "RELAY MUX , "RELAY ACTUATOR’ 

1080 REM 

1090 CALL SET.FUNCTION(DMM.01, DC VOLTS) 

1 100 CALL SET.RANGE(DMM.01,AUTOM) 

1110 CALL DIS ABLE.INT.TRIGGER(DMM.0 1 ) 

1 120 CALL ENABLE.OUTPUT(RELAY.MUX.01) 

1 130 CALL ENABLE.OUTPUT(RELAY.ACT.01) 

1 140 REM ************* PROGRAM TRAVERSE ******************* 

1 150 REM 

1 160 REM OPEN THE COM PORT AND INITIALIZE THE MOTOR SETTINGS 

1170 OPEN "COM 1 : 1200N T ,8,1 ,RS,CS ,DS ,CD" AS #1 

1180 REM SET MOTOR DEFAULT VALUES 

1190 DATA 2000,2000,2000,2,2,2,0.000125,0.000125,0.000125 

1200 RE.AD V1,V2,V3,R1,R2,R3,C1,C2,C3 

1210 REM DEFINE CHARACTERS FOR DATA REDUCTION ALGORITHM 
1220 RN2S=”RENAME AtRAW.DAT " 

1230 HEAD IS = " # X Y PI P2 P3 P4 P5 YAW ” 

1240 FORMATS= "## ##.## ##.## ###.#*# ###.### ###.### ###.### ###.### ###.##" 



1250 PRINT 

1260 PRINT "* **********************************************” 

1270 PRINT "** USER MUST SELECT "CAPS LOCK' FUNCTION **” 

1280 PRINT ******************************************** **★**" 

1290 PEM DISPLAY MOTOR DEFAULT SETTINGS 

***********************************•' 

INITIALIZED VALUES FOR ALL MOTOR SETTINGS:’' 
VELOCITY = 1000 STEPS/SEC" 

RAMP(MOTOR ACCELERATION) = 2 (6000 STEPS/SEC A 2)" 

DEFAULT INCREMENTAL UNITS ARE INCHES" 
************************************ 



1300 PRINT ” 
1310 PRINT " 
1320 PRINT ” 
1330 PRINT " 
1340 PRINT " 
1350 PRINT " 



1360 PRINT 

1370 PRINT "NOTE! ! USE MANUAL CONTROL TO INITIALIZE PROBE POSITION BEFORE" 
1380 PRINT" SELECTING COMPUTER CONTROLLED MOVEMENT. ’’ 

1390 PRINT 

1400 INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER 'MAN' OR CP')";CONS 
1410 IF CON$="CP" THEN 3490 

1420 REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY OR ACCELERATION 
RAMP 
1430 PRINT 
1440 PRINT 

1450 PRINT ’’ DO YOU WANT TO CHANGE THE VELOCITY OR ACCELERATION RAMP" 

1460 PRINT " DEFAULT SETTINGS? (Y OR N)" 

1470 PRINT 

1480 PRINT "IF "NO', THIS PROGRAM WILL THEN LET YOU DEFINE THE” 

1490 PRINT "DISTANCE YOU WANT TO MOVE (IN INCHES). IF "YES’," 

1500 PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS FOR ANY 
MOTOR." 

1510 PRINT 
1520 PRINT 
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1530 PRINT 

1540 INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? (Y OR N)";DS 
1550 IF D$="Y" THEN 1590 
1560 IF DS="N" THEN 2220 
1570 REM 

1580 REM **** OPERATOR SELECTED MOTOR VARIABLES ***** 

1590 PRINT 
1600 PRINT 

1610 INPUT "WHICH DEFAULT VALUE? (ENTER 'lTOR VELOC OR FOR ACCEL RAMP)";L 
1620 ON L GOTO 1690,1930 

1630 PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? (Y OR N)” 

1640 INPUT VS 

1650 IF V$="Y” THEN 1690 

1660 PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION RAMP? (Y OR N)" 

1670 IF RS="Y” THEN 1990 

1680 IF RS="N" THEN 1450 

1690 PRINT 

1700 PRINT 

1710 INPUT "WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2, OR 3)”;J 
1720 ON J GOTO 1730,1830,1880 
1730 PRINT 
1740 PRINT 

1750 INPUT "ENTER DESIRED VELOCITY OF MOTOR #1";V1 
1760 PRINT 
1770 PRINT 
1780 PRINT 

1790 PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? (Y OR N)" 

1800 INPUT VS 

1810 IF VS="Y" THEN 1690 

1820 IF VS="N” THEN 1430 

1830 PRINT 

1840 PRINT 

1850 INPUT "ENTER DESIRED VELOCITY OF MOTOR 2";V2 

1860 PRINT 

1870 GOTO 1780 

1880 PRINT 

1890 PRINT 

1900 INPUT "ENTER DESIRED VELOCITY OF MOTOR #3";V3 

1910 PRINT 

1920 GOTO 1780 

1930 PRINT 

1940 PRINT 

1950 INPUT "WHICH MOTOR ACCEL RAMP IX) YOU WANT TO CHANGE? (1, 2, OR 3)";K 
1960 ON K GOTO 1970,2060,2120 
1970 PRINT 
1980 PRINT 

1990 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #1”;R1 
2000 PRINT 
2010 PRINT 

2020 PRINT "DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER MOTOR? (Y OR 
N)?" 

2030 INPUT RMS 
2040 IF RMS="Y" THEN 1930 
2050 IF RMS="N" THEN 1450 
2060 PRINT 
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2070 PRINT 

2080 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #2 ’;R2 

2090 PRINT 

2100 PRINT 

21 10 GOTO 2000 

2120 PRINT 

2130 PRINT 

2140 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #3";R3 

2150 PRINT 

2160 PRINT 

2170 GOTO 2000 

2180 REM 

2190 REM DEFINE DISTANCE TO MOVE MOTOR 
2200 PRINT 
2210 PRINT 
2220 PRINT 

2230 REM INITIALIZE MOTOR INCREMENTS TO ZERO 

224011=0 

225012=0 

2260 13=0 

2270 PRINT 

2280 PRINT ” ********************************************************************" 
2290 PRINT " ** DEFINE WHICH MOTOR YOU WANT TO MOVE **" 

2300 PRINT ” ** **” 

2310 PRINT " ** NOTE!!! A POSITIVE (V) INCREMENT TO A MOTOR **" 

2320 PRINT " ** MOVES TRAVERSER AWAY FROM THAT PARTICULAR MOTOR **" 
2330 PRINT ” ** **" 

2340 PRINT ” ** - MOTOR #1 MOVES THE PROBE UPSTREAM AGAINST THE FLOW **" 
2350 PRINT ” ** - MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS WINDOW **" 
2360 PRINT '• ** - MOTOR #3 MOVES THE PROBE VERTICALLY DOWNWARD **" 
2370 PRINT " ****************************************************************** " 
2380 PRLNT 
2390 PRINT 

2400 INPUT "WHICH MOTOR DO YOU WANT TO MOVE? (1,2, OR 3)’’;L 
2410 ON L GOTO 2420,2680,2970 
2420 PRINT 
2430 PRINT 

2440 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 

2450 PRLNT ” ********* (ENTER DISTANCE IN INCHES) *********" 

2460 INPUT II 
2470 PRINT 

2480 PRINT” *********************************” 

2490 PRINT 

2500 PRINT "SUMMARY OF OPERATOR INPUTS 

2510 PRINT" MOTOR #1 VELOCITY = "; VI 

2520 PRINT " ACCELERATION RAMP = ”;R 1 

2530 PRINT " INCREMENTAL DISTANCE = ";I1;”INCHES" 

2540 PRINT” *********************************" 

2550 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y OR N)" 

2560 PRINT 

2570 PRINT "ENTER N’ TO START MOTOR MOVEMENT. ENTER ’Y' TO RETURN" 

2580 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

2590 INPUT VS 

2600 IF V$="Y" THEN 1430 

2610 GOSUB 3410 
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2620 PRINT 

2630 PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y OR N)?" 

2640 INPUT CS 

2650 IF CS="Y” THEN 2220 

2660 IF CS="N" THEN 3260 

2670 PRINT 

2680 PRINT 

2690 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #2?” 

2700 PRINT " ********* (ENTER DISTANCE IN INCHES) *********" 

2710 INPUT 12 
2720 PRINT 
2730 PRINT 

2740 REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
2750 PRINT" *********************************" 

2760 PRINT 

2770 PRINT "SUMMARY OF OPERATOR INPUTS:" 

2780 PRINT " MOTOR #2 VELOCITY = ";V2 

2790 PRINT " ACCELERATION RAMP = ";R2 

2800 PRINT " INCREMENTAL DISTANCE = ";I2;"INCHES" 

2810 PRINT" *********************************" 

2820 PRINT 
2830 PRINT 

2840 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y OR N)" 
2850 PRINT 

2860 PRINT "ENTER N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN- 
2870 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

2880 INPUT V$ 

2890 IF VS="Y" THEN 1430 
2900 GOSUB 3410 
2910 PRINT 

2920 PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y OR N)?” 

2930 INPUT CS 

2940 IF CS=”Y" THEN 2220 

2950 IF CS=”N" THEN 3260 

2960 PRINT 

2970 PRINT 

2980 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #3?" 

2990 PRINT " ********* (ENTER DISTANCE IN INCHES) *********" 

3000 INPUT 13 
3010 PRINT 
3020 PRINT 

3030 REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
3040 PRINT" *********************************” 

3050 PRINT 

3060 PRINT "SUMMARY OF OPERATOR INPUTS:” 

3070 PRINT " MOTOR #3 VELOCITY = ";V3 

3080 PRINT " ACCELERATION RAMP = ";R3 

3090 PRINT ” INCREMENTAL DISTANCE = ";I3;"INCHES" 

3100 PRINT 

3110 PRINT" *********************************" 

3120 PRINT 
3130 PRINT 

3140 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y OR N)” 
3150 PRINT 

3160 PRINT "ENTER N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN" 
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3170 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

3180 INPUT VS 

3190 IF V$="Y" THEN 1430 

3200 GOSUB 3410 

3210 PRINT 

3220 PRINT 

3230 INPUT "DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT (Y OR 
N)";M$ 

3240 IF M$=”Y" THEN 2210 
3250 PRINT 

3260 PRINT "DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR MOVEMENT?" 

3270 PRINT " ********* NOTE!!' ********* " 

3280 PRINT " ALL PREVIOUS MOTOR INCREMENT LNPUTS HAVE BEEN ZEROIZED." 

3290 PRINT "PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED MOVEMENT." 
3300 PRINT "***** (IF NO’, THE PROGRAM WILL END). *****" 

3310 PRINT 

3320 INPUT "DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y OR N)";N$ 

3330 IF NS="Y" THEN 3500 

3340 PRINT 

3350 PRINT 

3360 PRINT 

3370 PRINT ” ****************************" 

3380 PRINT " THE PROGRAM HAS ENDED." 

3390 PRINT ” ****************************" 

3400 END 

3410 REM ******* MOTOR MOVEMENT SUBROUTINE ********* 

3420 PRINT # 1, "&" :PRINT #1, "E";"C1=";C1;":C2=";C2;":C3=";C3 
3430 PRINT #1, "1 1=";1 1;":V1=”;V1;":R1=";R 1 ; 

3440 PRINT #1, ":I2=";I2;":V2=";V2;":R2=":R2 
3450 PRINT #1, "I3=";I3;":V3=";V3;":R3=";R3;”:@" 

3460 RETURN 

3470 REM ************************ ******************** 

3480 REM **** a t ************************** ************* 

3490 PRINT 

3500 REM ******* COMPUTER CONTROLLED MOVEMENT ******* 

3510 PRINT 

3520 PRINT "THE PRESSURE DATA WILL BE WRITTEN TO FILES ON DRIVE ’A' " 

3530 PRINT 

3540 PRINT "YOU WILL BE ASKED TO INPUT FILE NAMES FOR THESE." 

3550 PRINT 

3560 INPUT "IS A FORMATTED DISK IN DRIVE 'A'? PRESS 'ENTER' TO CONTINUE";D$ 

3570 PRINT 
3580 PRINT 
3590 PRINT 
3600 PRINT' 

3610 PRINT 
3620 PRINT ’ 

3630 PRINT ’ 

3640 PRINT ' 

3650 PRINT ’ 

3660 PRINT ' 

3670 PRINT ' 

3680 PRINT 

3690 REM SET INITIAL MOVEMENT DISTANCE AND NUMBER OF DATA POINTS TO ZERO 
3700 HT=0 



******************* *•**•••*•••*•****•******••**•< 

** NOTE !!! **" 

* * COMPUTER CONTROLLED MOVEMENT * * ” 

** IS PROGRAMMED WITH A **" 

** DEFAULTED NEGATIVE MOTOR INCREMENT **" 
** (LE. MOTOR #3 WILL MOVE UPWARD **" 

** BY ENTERING A (+) DISTANCE). 
************************************************" 
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3710 WD=0 
3720 DIST=0 
3730 XPT=0 
3740 YPT=0 
3750 N=0 
3760 PRINT 
3770 PRINT 

3780 INPUT "WHAT IS THE DIMENSION ( X , Y ) (IN INCHES) THAT YOU WANT TO 
MEASURE." ;WD,HT 
3790 PRINT 

3800 INPUT "WHAT IS THE STEP (IN INCHES) THAT YOU WANT TO MOVE.";DIST 

3810 YPT=INT(HT /DIST) + 1 

3820 XPT=INT(WD /DIST)+ 1 

3830 N=XPT*YPT 

3840 PRINT 

3850 PRINT "THERE ARE ";XPT;" * ";YPT;" = ";N;" POINTS TO BE MEASURED " 

3860 PRINT 

3870 INPUT "ARE THE NUMBER OF POINTS IS OK.(Y OR N)";C$ 

3880 IF CS="N” THEN 3780 
3890 CLS 
3900 N=XPT 

3910 IF (N < 1) OR (N > 99) GOTO 3780 

3920 REM *** GENERATING STRING STRING SEGMENTS FOR DATA FILE NAMES 
3930 BS = MIDS(STRS(1), 2): REM ** STRING NUMBER ”1" 

3940 ES = MIDS(STRS(N), 2): REM ** ENDING STRING NUMBER "N” 

3950 XS = "XXXXXX" 

3960 EXS = ".DAT” 

3970 CLS 

3980 PRINT "DATA FILES WILL BE INCREMENTED FROM:" 

3990 PRINT 

4000 PRINT (XS + BS + EXS); " TO "; (XS + ES + EXS) 

4010 PRINT 
4020 PRINT 

4030 INPUT "ENTER DATA FILE NAME (6 CHARACTERS MAX - NO EXTENSION)”;F2S 
4040 PRINT 
4050 PRINT 

4060 IF LEN(F2S) > 6 OR LEN(F2S) < 1 GOTO 4030 
4070 CLS 

4080 PRINT N; "DATA FILES WILL BE GENERATED AND INCREMENTED AS FOLLOWS 
4090 PRINT 
4100 PRINT 

4110 PRINT (F2S + BS + EXS); " TO ”; (F2S + ES + EXS) 

4120 PRINT 
4130 PRINT 

4140 INPUT "ARE THE NUMBER OF POINTS AND FILE NAMES OK.(Y OR N)”; CS 

4150 IF CS = "N" GOTO 3780 

4160 IF CS = "Y" GOTO 4180 

4170 GOTO 4140 

4180 CLS 

4190 PRINT 

4200 PRINT 

4210 REM SET INITIAL POSITION DATA 

4220 X(1)=-DIST 

4230 Y(1)=-DIST 

4240 FOR IX=2 TO XPT+1 
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4250 X(IX)=0 

4260 NEXT IX 

4270 FOR JY=2 TO YPT+1 

4280 Y(JY)=0 

4290 NEXT JY 

4300 FOR 1=1 TO XPT 

430211=0 

4304 12=0 

4306 13=0 

4310 FOR J=1 TO YPT 

4320 REM MOTOR CP-CONTROLLED MOTOR MOVEMENT 

433011=0 

434012=0 

4350 13=0 

4360 REM EACH POINT TAKE 10 TIMES READINGS 
4370 X(I+ 1 )=X(I)+DIST 
4380 XPT(J)=X(I+1) 

4390 Y(J+ 1 )= Y (J)+DIST 
4400 YPT(J)=Y(J+1) 

4405 INPUT " ADJUST THE WHEEL TO MAKE THE P2 =P3, INPUT THE YAW ANGLE";YAW(J) 
4408 PRINT 

4410 INPUT " PRESS 'ENTER' TO ST/ RT THE MEASUREMENT" ;MO VES 
4420 REM 

4430 REM READ FIVE CHANNELS AND DISPLAY THE DATA 

4440 REM 

4450 STEPPERS 

4460 SWITCH = 3 

4470 HOMER=8 

4480 DELAY 1 = .1 

4490 DELAY2 = 1 

4500 REM SET THE S.V PORT TO #4 
4510 FOR IL=1 TO 3 
4520 THYME = TIMER 

4530 CALL OUTPUT(RELAY. ACT.01, STEPPER) 

4540 CHKTIME = TIMER 

4550 IF CHKTIME < (THYME + DELAY1) GOTO 4540 
4560 CALL OPEN.CHANNEL(RELAY. ACT.01, SWITCH1 
4570 CLS 
4580 NEXT IL 
4590 PRINT 

4600 PRINT " NOW IS ";J;” POINT ” 

4610 REM START MEASURE FROM PORT 4 TO PORT 8 
4620 FOR JJ=1 TO 5 

4630 CALL OUTPUT(REL AY. ACT.01, STEPPER) 

4640 CHKTIME = TIMER 

4650 IF CHKTIME < (THYME + DELAY 2) GOTO 4640 
4660 REM EACH PORT SAMPLE 10 TIMES 
4670 FOR 11=1 TO 10 
4680 ROUT= 1 

4690 CALL OUTPUT(RELAY.MUX.01,ROUT) 

4700 CALL ME ASURE(DMM.01 .VOLTS) 

4710 PA(DJJ)= VOLTS 
4720 NEXT II 

4730 CALL OPEN.CHANNEL(RELAY. ACT.01 .SWITCH) . 

4740 IF JJ=5 THEN 4760 
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4750 NEXT JJ 

4760 REM HOME THE S.V PORT TO #48 
4770 CALL OUTPUT(REL AY. ACT.01, HOMER) 

4780 CALL OPEN.CHANNEL(RELAY. ACT.01, HOMER) 

4790 REM 

4800 REM DISPLAY THE SAMPLE DATA 
4810 REM 

4820 PRINT HEAD1S 
4830 FOR IS= 1 TO 10 

4840 PRINT USING FORMATS;IS,XPT(J),YPT(.r)PA(IS,l)PA(IS.2)PA(IS,3)PA(IS,4), 
PA(IS,5),Y AW(J) 

4850 NEXT IS 
4860 REM 

4870 REM AVERAGE THE DATA 
4880 REM 

4890 FOR JA = 1 TO 5 

4900 TOTAL = 0 

4910 FOR IA = 1 TO 10 

4920 TOTAL = TOTAL + PA(IA JA) 

4930 NEXT IA 

4940 AVERAGE = TOTAL /10 
4950 P(JA)= AVERAGE 
4960 NEXT JA 
4970 PRINT 

4980 PRINT "THE AVERAGE ARE: ” 

4990 PRINT 
5000 PRINT HE AD IS 
5010 FOR JD=1 TO 5 
5020 PP(J JD)=P(JD) 

5030 NEXT JD 

5040 PRINT USING FORMATS J ,XPT ( J ) , YPT ( J) ,PP( J , 1 ) PP(J ,2) PP( J ,3) PP( J ,4) ,PP(J ,5) , Y A W(J) 
5045 PRINT 

5050 PRINT "DO YOU WANT RE-MEASURE AGAIN (Y / N)" 

5060 PRINT 

5062 PRINT "IF Y' WILL RE-SAMPLE AGAIN." 

5064 PRINT 

5070 INPUT "IF N’ WILL MOVE THE TRAVERSER STEP UPWARD (WAIT 7 SEC )";CS 
5075 PRINT 

5080 IF CS="Y" THEN 4405 
5082 IF CS="N" THEN 5090 
5084 GO TO 5070 
5090 IF J=YPT THEN 5160 
5100 REM 

51 10 REM MOVE THE TRAVERSER STEP UPWARD. 

5120 REM 
5130 13=-DIST 
5140 GOSUB 3410 
5150 NEXT J 

5160 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

5170 OPEN "A:\RAW.DAT" FOR OUTPUT AS #2 
5180 PRINT #2 .HE ADIS 
5190 FOR ID=1 TO YPT 

5200 PRINT #2 .USING FORMATS;ID,XPT(ID),YPT(ID)PP(ID,l)PP(ID,2)PP(ID,3). PP(ID,4), 
FP(ID,5),YAW(ID) 

5210 NEXT ID 
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5220 CLOSE #2 

5230 REM *** GENERATING INCREMENTED DATA FILE NAME 
5240 IF (I > 10) OR (I = 10) THEN IS = MID$(STR$(I), 2) 

5250 IF (I < 10) THEN IS = (MIDS(STRS(0), 2) + MIDS(STR$(I), 2)) 

5260 FI2S = (F2S + IS + EXS) 

5270 PRINT 

5280 PRINT " WRITING DATA FILE FI2S 

5290 DF2S=RN2S+FI2S 

5300 REM ** RENAME DATA FILE 

5310 SHELL DF2S 

5320 REM 

5330 REM MOVE THE TRAVERSER TO THE NEXT SAMPLE POSITION 
5340 REM 
5350 PRINT 

5360 IF I=XPT THEN 5430 

5370 INPUT "THEN PRESS ’ENTER' FOR NEXT COLUMN SAMPLE( 90 SEC) ";MOVES 

5390 I2=-DIST 

5400 I3=HT 

5410 GOSUB 3410 

5420 NEXT I 

5430 CLS 

5440 PRINT "ALL MOVEMENTS COMPLETE" 

5450 PRINT 
5460 PRINT 

5470 PRINT "YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER MOVEMENT ( Y OR N )?" 
5480 PRINT 

5490 PRINT "IF ’Y\ THE PROGRAM WILL TAKE YOU TO MANUAL CONTROL SUBROUTINE." 
5500 PRINT "IF K', THE PROGRAM WILL END." 

5510 PRINT 

5520 INPUT "ANOTHER MOVEMENT" ;RS 
5530 IF RS = "Y" THEN 1370 
5540 IF RS = "N" THEN 3370 
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APPENDIX B. CALP PROGRAM 



1 DEF SEG.CLEAR ,&HFEOO:GOTO 4 'BEGIN PCIB PROGRAM SHELL 

2 GOTO 1000 * USER PROGRAM 

3 GOTO 900 ' ERROR HANDLING 

4 I=&HFE00 ’ COPYRIGHT HEWLETT-PACKARD 1984,1985 

5 PCIB .DIRS=ENVIRON$(’ , PCIB ") 

6 IS=PCIB.DIRS+"\PCIBILC.BLD" 

7 BLOAD IS,I 

8 CALL I(PCIB.DIRS4% J%):PCIB.SEG=I% 

9 IF J%=0 THEN GOTO 13 

10 PRINT "UNABLE TO LOAD.’’; 

1 1 PRINT ” (ERROR 

12 END 

13 ’ 

14 DEF SEG=PCIB.SEG:O.S=5:C.S=10:I.V=15 

15 I.C=20:L.P=25:LDJTLE=30 

16 GET.MEM=35:L.S=40:PANELS=45:DEF.ERR=50 

17 PCIB.ERRS=STRINGS(64,32) : PCIB.NAMES=STRING$(16,32) 

18 CALL DEF.ERR (PCIB. ERR .PCIB .ERRS,PCIBNAMES,PCIB.GLB ERR) : PCIB.BASERR=255 

19 ON ERROR GOTO 3 

20 J=-l 

21 IS=PCIB.DIRS+"\PCIB.SYN" 

22 CALL O.S(IS) 

23 IF PCIB. ERR <>0 THEN ERROR PCIB.BASERR 
241=0 

25 CALL I.V(IJIEAD.REGISTER,READ.SELFID,DEFTNE,INITIALIZE.SYSTEM) 

26 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

27 CALL I.V(I,ENABLE.SYSTEM,DISABLE.SYSTEM,INrnALIZE,POWER.ON) 

28 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

29 CALL LV(I,MEASURE,OUTPUT,START,HALT) 

30 IF PCIB. ERR <>0 THEN ERROR PCIB.BASERR 

3LCALLLV(I.ENABLE.INT.TRIGGER,DISABLE.INT.TRIGGER,ENABLE.OUTPUT,DISABLE.OUTT) 

32 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

33 CALL I. VG,CHECK.DONE,GET.STATUS, SET JUNCTION .SET.RANGE) 

34 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

35 CALL I.V(I,SET.MODE,WRITE.CAL JEAD.CAL.STORE.CAL) 

36 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

37 CALL I.V(I .DELAY, SAVE.SYSTEMJJ) 

38 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

39 1=1 

40 CALL I.V(I,SET.G ATE TI ME, SET.SAMPLES,SET. SLOPE, SET. SOURCE) 

41 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

42 CALL I.C(I FREQUENCY, AUTO.FREQ JERIOD.AUTO.PER) 

43 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

44 CALL I.C(I .INTERVAL .RATIO, TOTALIZER 100 MILL I) 

45 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

46 CALL I.C(IJIJIOJIOOJIKILO) 

47 IF PCIB JRRoO THEN ERROR PCIB.BASERR 

48 CALL I.CG J IOMEGA J100MEGA,CHAN.A,CHAN.B) 

49 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

50 CALL I.Ca, POSITIVE, NEGATIVE,COMN,SEPARATE) 
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51 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

521=2 

531=3 

54 CALL I.V(I,ZERO.OHMS,SET.SPEED,J,J) 

55 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

56 CALL I.C(1 ,DC V OLTS ,AC V OLTS .OHMS ,R200MILLI) 

57 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

58 CALL I.C(I,R2,R20,R200,R2iaLO) 

59 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

60 CALL I.C(I,R20KILO,R200KILO,R2MEGA,R20MEGA) 

61 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

62 CALL I.C(l,AUTOM,R2.5,R12.5,J) 

63 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 
641=4 

65 CALL I.V(I,SET.COMPLEMENT,SET.DRIVER, OUTPUT.NO. WAIT.ENABLE.HANDSHAKE) 

66 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

67 CALL I.V(1,DISABLE.HANDSHAKE,SET.THRESH0LD,SET.START.BIT,SET.NUM.BITS) 

68 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

69 CALL I.VCLSET.LOGIC.SENSEJJJ) 

70 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

71 CALL I.C(I .POSITIVE .NEGATIVE, TWOS, UNSIGNED) 

72 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

73 CALL I.C(I,OC,TTL,RO,R1) 

74 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

75 CALL I.C(I,R2,R3,R4,R5) 

76 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

77 CALL I.C(I,R6,R7,R8,R9) 

78 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

79 CALL I.C(I,R10,R1 1,R12,R13) 

80 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

81 CALL I.C(I,R14,R15,R16,J) 

82 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

83 1=6 

84 CALL I.V(I,SET.FREQUENCY,SET.AMPLITUDE,SET.OFFSET,SET.SYMMETRY) 

85 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

86 CALL I.VG.SET.BURST.COUNTJJJ) 

87 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

88 CALL I.C(I,SINE, SQUARE, TRIANGLE.CONTINUOUS) 

89 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

90 CALL I.C(I,GATED,BURSTT J) 

91 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
921=7 

93 CALL I.V(I AUTOSCALE, CALIBRATE, SET.SENSmVITY.SET.VERT.OFFSET) 

94 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

95 CALL I.V(I,SET.COUPLING,SET.POLARITY,SET.SWEEPSPEED,SET.DELAY) 

96 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

97 CALL I. V(I, SET .TRIG .SOURCE.SET .TRIG. SLOPE.SET.TRIG.LEVEL.SET.TRIG. MODE) 

98 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

99 CALL I.V(I,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO,GET.TIMEBASE.INFO) 

100 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

101 CALL I.V(I,GET.TRIG.INFO,CALC.WFVOLT,CALC.WFTIME,CALC.WF.STATS) 

102 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

103 CALL I. V(1,CALC.RISETIME,CALC.FALLTIME,CALC. PERIOD .CALC.FREQUENCY) 

104 IF PCIB.ERR 0 O THEN ERROR PCIB.BASERR 

105 CALLI.V(I,CALC.PLUSWIDTH,CALC.MINUSWIDTH,CALC.OVERSHOOT,CALC.PRESHOOT) 
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106 IF PCIBERRoO THEN ERROR PCIB.BASERR 

107 CALL I. V(I,CALCPK.TO.PK,SET.TIMEOUT .SCOPE. START .MEASURE. SINGLE. WF) 

108 IF PCIBERRoO THEN ERROR PCIB.BASERR 

109 CALL I.VG.MEASURE.TWO.WFJJ J) 

110 IF PCIBERRoO THEN ERROR PCIB.BASERR 

111 CALL I.COPlONANOElOONANOElMICROjUOMICRO) 

112 IF PCIBERRoO THEN ERROR PCIB.BASERR 

113 CALL I ,C(I E 1 OOMICRO E 1 MILLI ,R 10MILLI E 1 00MILLI) 

114 IF PCIBERRoO THEN ERROR PCIB.BASERR 

115 CALL I.CaElE10E20NANOE200NANO) 

116 IF PCIBERRoO THEN ERROR PCIB.BASERR 

117 CALL I.CGE2MICROE20MICROE200MICROE2MILLI) 

118 IF PCIBERRoO THEN ERROR PCIB.BASERR 

119 CALL I.C(1E20MILLIE200MILLIE2E20) 

120 IF PCIBERRoO THEN ERROR PCIB.BASERR 

121 CALL I.C(IE50NANOE500NANOE5MICROE50MICRO) 

122 IF PCIBERRoO THEN ERROR PCIB.BASERR 

123 CALL I.CGE500MICROE5MILLIE50MILLIE500MILLI) 

124 IF PCIBERRoO THEN ERROR PCIB.BASERR 

125 CALL I.C(IE5E50,CHAN.A,CHAN.B) 

126 IF PCIBERRoO THEN ERROR PCIB.BASERR 

127 CALL I.C(I£XTERNALPOSITI\E,NEGATIVE,AC) 

128 IF PCIBERRoO THEN ERROR PCIB.BASERR 

129 CALL LC0JX\TRIGGERED r AUTO.TRIG r AUTOEEVEL) 

130 IF PCIBERRoO THEN ERROR PCIB.BASERR 

131 CALL I .C(I,X 1 ,X 1 0,ST ANDARD, AVER AGE) 

132 IF PCIBERRoO THEN ERROR PCIB.BASERR 

133 1=8 

134 CALL I. V(I,OPEN.CHANNEL, CLOSE. CHANNEL JJ) 

135 IF PCIBERRoO THEN ERROR PCIB.BASERR 

136 CALL C.S 

137 IF PCIBERRoO THEN ERROR PCIB.BASERR 

138 IS=PCIB.DIR$+"\PCIB.PLD" 

139 CALL L.P(IS) 

140 IF PCIBERRoO THEN ERROR PCIB.BASERR 

141 IS="DMM.01":I=3:J=0:K=0:L=1 

142 CALL DEFTNE(DMM.01,IS,U,KE) 

143 IF PCIBERRoO THEN ERROR PCIB.BASERR 

144 IS="FUNC.GEN.0 1 ”:I=6: J=0:K= 1 : L= 1 

145 CALL DEFINE(FUNC.GEN.01 JS.U.KE) 

146 IF PCIBERRoO THEN ERROR PCIB.BASERR 

147 IS="SCOPE.01":I=7:J=0:K=2:L=1 

148 CALL DEFINE(SCOPE.01,I$J r F,K,L) 

149 IF PCIBERRoO THEN ERROR PCIB.BASERR 

150 IS="COUNTER.01":I=1:J=0:K=3:L=1 

151 CALL DEFINE(COUNTER.01,ISJJ,KE) 

152 IF PCIBERRoO THEN ERROR PCIB.BASERR 

153 I$="DIG.IN.0r , :I=4:J=0:K=4:L=l 

154 CALL DEFINE (DIG. IN.01, IS, I.J.K.L) 

155 IF PCIBERRoO THEN ERROR PCIB.BASERR 

156 I$="DIG.OUT.01":I=4:J=1:K=4:L=1 

157 CALL DEFINE(DIG.OUT.01 JS.IJ.KL-) 

158 IF PCIBERRoO THEN ERROR PCIB.BASERR 

159 IS="RELAY.ACT.01":I=8:J=0:K=5:L=1 

160 CALL DEFTNE(RELA Y.ACT.0 1 JS,I,J,KE) 
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161 IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

162 IS="RELAY.MUX.01":I=2:J=0:K=6:L=1 

163 CALL DEFTNE(RELAY.MUX.01,I$4,J,K,L) 

164 IF PCIBPRRoO THEN ERROR PCIB.BASERR 

800 IS=ENVIRON$("PANELS")+"\PANELSJEXE" 

801 CALL L.S(I$) 

899 GOTO 2 

900 IF ERR=PCIB .B ASERR THEN GOTO 903 

901 PRINT "BASIC ERROR #”;ERR;" OCCURRED IN LINE ";ERL 

902 STOP 

903 TMPERR=PCIB.ERR:IF TMPERR=0 THEN TMPERR=PCIB.GLBERR 

904 PRINT "PC INSTRUMENT ERROR #”;TMPERR;" DETECTED AT LINE ";ERL 

905 PRINT "ERROR: ";PCIB.ERRS 

906 IF LEFTS (PCIB .NAMES, 1)<>CHRS(32) THEN PRINT "INSTRUMENT: ";PCIB.NAME$ 

907 STOP 

908 COMMON PCIB.DIRS.PCIB.SEG 

909 COMMON LD.FILE,GET.MEMPANELS,DEF.ERR 

910 COMMON PCIB.BASERR,PCIB.ERRJPCIB.ERRS,PCIBJn , AMES,PCIB.GLBERR 

91 1 COMMON READ JtEGISTER,READ.SELFID, DEFINE, INITIALIZE.SYSTEM, 

ENABLE.SYSTEM,DISABLE.SYSTEM,INITIALIZE,POWER.ON,MEASURE,OUTPUT, 

START.HALTJENABLE.INT.TRIGGER.DISABLE.INT.TRIGGER^NABLE.OLTPUT, 

DISABLE.OUTPUT,CHECK.DONE,GET.STATUS 

912 COMMON SET.FUNCTION,SET.RANGE,SET.MODE,WRITE.CAL,READ.CAL, 

STORE.CAL,DELAY,SAVE.SYSTEM,SET.GATETIME,SET.SAMPLES,SET.SLOPE, 
SET.SOURCE,ZERO.OHMS,SET.SPEED,SET.COMPLEMENT,SET.DRIVER, 
OUTPUT.NO. WAITLNABLE.HANDSHAKE, DISABLE. HANDSHAKE 

913 COMMON SET.THRESHOLD,SET.START.BIT,SET.NUM.BITS,SET.LOGlC. SENSE, 

SET.FREQUENCY,SET.AMPLITUDE,SET.OFFSET,SET.SYMMETRY, 
SET.BURST.COUNT, AUTOSCALE, CALIBRATE, SET.SENSIT1VITY, 
SET.VERT.OFFSET,SET.COUPLING,SET.POLARITY,SET.SWEEPSPEED 

914 COMMON SET.DELAY.SET.TRIG. SOUR CE.SET. TRIG. SLOPE, SET. TRIG.LEVEL, 

SET.TRIG.MODE,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO, 

GET.TlMEBASE.INFO,GET.TRIG.INFO,CALC.WFVOLT,CALC.WFTIME, 

CALC.WF.STATS,CALC.RISETIME,CALC.FALLTIME,CALC.PERIOD 

915 COMMON CALC.FREQUENCY,CALC.PLUSWIDTH,CALC.MINUS WIDTH, 

C ALC.O VERS HOOT.CALC PRES HOOT, CALC.PK.TO.PK.SET.TIMEOUT, 

SCOPE.START,MEASURE.SINGLE.WF,MEASURE.TWO.WF,OPEN.CHANNEL, 

CLOSE.CHANNEL 

916 COMMON FREQUENCY, AUTO.FREQPERIOD,AUTOPER,INTERVAL,RAT10, 

TOTALIZEP 1 00MILLI.R 1 P 1 0,R 1 OOP 1 KILO.R 1 0MEG AP 1 00MEG A,CHAN. A, 
CHAN.BPOSITIVE,NEGATIVE,COMN, SEP AR ATE, DC VOLTS ,AC VOLTS, OHMS, 
R200MILLIP2P20P200P2KILOP20KILOP200KILO 

917 COMMON R2MEGAP20MEGA,AUTOMP2.5P12.5POSITIVE,NEGATIVE,TWOS, 

UNSIGNED,OC,TTLPO,R1P2,R3P4,R5P6P7P8P9P10P11P12P13,R14,R15, 

R1 6,SINE,SQU ARE.TRIANGLE, CONTINUOUS ,GATED,BURST P 10NANO, 
RIOONANOPIMICROPIOMICROPIOOMICRO 

918 COMMON R1MILLIP10MILLIP100MILLIP1P10P20NANO,R200NANO,R2MICRO, 

R20MICROP200MICROP2MILLIP20M]LLIP200MILLI,R2,R20P50NANO, 
R500NANOP5MICROP50MICROP500MICROP5MILLIP50MILLIP500MILLIP5, 
R50.CHAN. A.CHAN.B EXTERNAL POSITIVE 

919 COMMON NEGATIVE,AC,DC,TRIGGERED,AUTO.TRIG,AUTO.LEVEL,X1,X10, 

ST AND ARD AVER AGE 

920 COMMON DMM.01pUNC.GEN.01,SCOPE.01,COUNTER.01,DIG.IN.01,DIG.OUT.01, 
RELAY. ACT.01 PELAY.MUX.01 

999 LND PCIB PROGRAM SHELL 
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1000 REM THIS STEP INTTIALZES THE HP SYSTEM 
1010 REM THIS STEP INTTIALZES THE HP SYSTEM 
1020 CLS 

1030 OPTION BASE 1 

1040 DIM P(10) J > A(50,6),PP(50,6) V XPT(40),CAL(40) 

1050 CALL INITIALIZED YSTEM(PGMSHEL.HPC) 

1060 REM 

1070 REM ALL PC DEVICES NOW HAVE AN INITIAL STATE 
1080 REM SET FUNCTION ON THE DMM AND RELAY MUX 
1090 REM 

1 100 CALL SET.FUNCTIONCDMM.OIJOC VOLTS) 

1 1 10 CALL SET.RANGE(DMM.01,AUTOM) 

1120 CALL DISABLE.INT.TRIGGER(DMM.01) 

1 130 CALL ENABLE. OUTPUT(RELAY.MUX. 01) 

1 140 FflRMA’K-"## MM MMMM MM MMMM MM MMMM MM MMMM MM MM MM MM MMMM" 

1 I VylVlfl/i 1 ft 77 Tt TT . TT tttttT 77 Tt ,Tt7TTtTt TtTT .77 FT TT 77 TtTT « FI fl Tf Tj 71 n .TTTTrTTT TtTT 9 ftftftft 

1200 FOR 1=1 TO 10 
1210 CAL(I)=0.0 
1220 NEXT I 
1510 REM 

1520 REM READ THE VOLTAGE OF 48TH CHANNEL AND DISPLAY THE DATA 
1530 REM 

1540 PRINT " CHOOSE 6 POINTS" 

1550 PRINT 

1550 PRINT "THE CALIBRATION WILL BE STORES IN ’CAL.DAT" 

1560 REM 

1570 REM BEGIN SAMPLING LOOP 

1580 REM 

1600 FOR J=1 TO 1 

1610 PRINT 

1630 FOR JJ=1 TO 6 

1631 INPUT "INPUT THE CALIBRATION PRESSURE";CAL(JJ) 

1632 INPUT "PRESS 'ENTER' TO START MEASURE MENT";MOVES 
1640 FOR 11=1 TO 10 

1650 ROUT=l 

1660 CALL OUTPUT (RELAY. MUX. 01 .ROUT) 

1670 CALL ME ASURE(DMM.01, VOLTS) 

1680 PA(H JJ)=VOLTS 

1690 NEXT II 

1700 IF JJ=6 THEN 1740 

1730 NEXT JJ 

1740 REM 

1750 REM DISPLAY THE SAMPLE DATA 
1760 REM 

1780 FOR IS= 1 TO 10 

1790 PRINT USING FORMATS;IS,PA(IS,l),PA(IS,2),PA(IS,3),PA(IS,4),PA(IS,5),PA(IS,6) 
1800 NEXT IS 
1810 REM 

1820 REM AVERAGE THE DATA 
1830 REM 

1840 FOR JA = 1 TO 6 

1850 TOTAL = 0 

1860 FOR IA= 1 TO 10 

1870 TOTAL = TOTAL + PAG A J A) 

1880 NEXT I A 

1890 AVERAGE = TOTAL /10 
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1900 P(JA)= AVERAGE 
1920 NEXT J A 
1930 PRINT 

1940 PRINT "THE AVERAGE ARE: " 

2000 FOR JD=1 TO 6 
2010 PP(J JD)=P(JD) 

2020 NEXT JD 

2055 PRINT USING FORMAT$;JJ > P(J,l),PP(J^),PP(J,3)PP(J,4),PP(J,5) > PP(J,6) 
2070 PRINT 

2080 INPUT "DO YOU WANT RE-MEASURE AGAIN ? (Y / N)";CS 
2090 IF C$="Y" THEN 1580 

2101 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

2102 OPEN "A:\CAL.DAT" FOR OUTPUT AS #2 

2106 FOR ID=1 TO 6 

2107 PRINT #2,USING FORMATS;IDJ > P(J,ID),CAL(ID) 

2108 NEXT ID 

2109 CLOSE #2 
2210 NEXT J 
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APPENDIX C. CONVERT PROGRAM 



************************************************************************ 

THIS PROGRAM CONVERTS THE VOLTAGE OF TRANSDUCER INTO PHYSICAL 
PRESSURE, VELOCITY ,Y AW ANGLEJ’ITCH ANGLE.TOTAL PRES SURE .TOTAL PRESSURE 
COEFFICIENT, STATIC PRESSURE AND STATIC PRESSURE COEFFICIENT.THOSE DATA ARE 

USED TO PLOT VELOCITY MAPPING AND PRESSURE CONTOUR. 
************************************************************************ 

CHARACTER* 12 FNAME 
CHARACTER* 12 NAME 
CHARACTER*2 A(13) 

CHARACTER*80 ST 

REAL LO,PI,PF,TI,TF,K,E 

DATA A/’0r, , 02', , 03','O4 , ,’05','06','07’,’08','09', 

* T0',Tr,T2',T3'/ 

NAME='R0B 1A100.DAT' 

* CONVERT THE PRESSURE UNIT FROM INHG TO PSF 

PI=29.81 

PF=29.89 

PATM=(PI+PF)*2.54*27.845/2. 

K=0.889 
E=0.0123 
R= 1718 
TI=62 
TF=68 

T=(TI+TF)/2.+460 

LO=PATM/(R*T) 

* OPEN A NEW FILE TO STROE THE REDUCED DATA 

OPEN(2TILE='RESULTO.DAT,STATUS='NEW') 

* OPEN A SEQUENTIAL OF DATA FILE 

DO 20 1=1,13 
NAME(7:8)=A(I) 

FN AME=N AME 
OPEN( 1 f ILE=FNAME) 

READ(1,100)ST 
100 FORMAT(A65) 

15 READ(1,1000JEND=30)NO,X,Y,V1,V2,V3,V4,V5,BETA 

1000 FORMAT(I2,F7.2,F6.2,5F9.3,F8.2) 

* CONVERT THE VOLTAGE TO PRESSURE IN LBF/FT**2 

P1=DELTAP(V1)*2.0475+PATM 
P2=DELTAP(V2)*2.0475+PATM 
P3=DELT AP(V 3) ‘2.0475+PATM 
P4=DELT AP(V4)*2.0475+PATM 
P5=DELTAP(V5)*2.0475+PATM 

* CALCULATE THE PITCH ANGLE IN DEGREE 

P=(P4-P5)/(P1 -P2) 

ALPHA=FPITCH(P) 

* CALCULATE THE VELOCITY IN FT/SEC 

YSLOP=FY SLOP( ALPHA) 

VELM= S QRT ((2* YS LOP* (P 1 -P2))/(LO* K)) 

VEL=VELM*(1+E) 

* CALCULATE THE FREE STREAM AND LOCAL DYNAMIC PRESSURE 
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QM1=7.2*2.0475/K 

QM=LO*VEL**2/2. 

Q1=QM1*(1+2*E) 

Q=QM*(1+2*E) 

* CALCULATE THE YAW ANGLE IN DEGREE 

YAW=FYAW(BETA) 

* CALCULATE THE LOCAL TOTAL PRESSURE COEFFICIENT 

PTC=FPT( ALPHA) 

PT1=P1-Q*PTC 

PT=PT1/144 

CPT=(PT1 -PATM-Q 1 )/Q 1 

* CALCULATE THE LOCAL STATIC PRESSURE COEFFICIENT 

PS1=PT1-Q 
PS=PS1/144 
CPS=(PS 1-PATM)/Q1 

WRITE(2,2000)NO,X,Y,\ t EL,YAW,ALPHA,PT,CPTJ > S,CPS 
2000 FORMAT (15 ,9F1 0.3) 

GO TO 15 
30 CLOSE(l) 

20 CONTINUE 
CLOSE(2) 

STOP 

END 

*********** ************************************************** 

* THIS FUNCTION CONVERT THE VOLTAGE TO PHYSICAL PRESSURE 

FUNCTION DELTAP(X) 

REAL DELTAP.X 

DELTAP=0.0214+4.2321*X 

END 

************************************************************* 

* THIS FUNCTION CALCULATES THE PITCH ANGLE 

FUNCTION FPITCH(X) 

REAL FPITCH.X 

FPITCH=3.759+53.7568*X-1.3085*X**2-1.6583*X**3 
* -0.8061 *X**4+16. 51 15*X**5 

END 

************************************************************ 

* THIS FUNCTION CALCULATES THE VELOCITY PRESSURE COEFFICENT 

FUNCTION FYSLOP(X) 

REAL FYSLOP.X 
IF(X.LT.-10)THEN 

FYSLOP=0.981-0.0102*X-3.000E-4*X**2-2.500E-6*X**3 
ELSE IF((X.GE.-10).AND.(XUE.10))THEN 
FYSLOP=0.98-0.006*X+2.000E-4*X**2 
ELSE 

FYSLOP=0.9801-0.0035*X-1.143E-4*X**2+5.833E-6*X**3 
END IF 
END 

************************************************************* 

* THIS FUNCTION CALCULATES THE YAW ANGLE 

FUNCTION FYAW(X) 

REAL FYAW.X 

IF((X.GE.0).AND.(X.LE.180)) THEN 
FYAW=-X 
ELSE 
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FYAW=360-X 
END IF 
END 

************************************************************* 

* THIS FUNCTION CALCULATES THE TOTAL PRESSURE COEFFICIENT 
FUNCTION FPT(X) 

REAL FPT,X 
EF(XLE.-30) THEN 
FPT=-0.01 

ELSE IF((X.GT.-30).AND.(X.LT.-20)) THEN 
FPT=0.02+1.00E-3*X 

ELSE IF((X.GE.-20).AND.(XEE.30)) THEN 
FPT=0 
ELSE 

FPT=0.03-1.00E-3*X 
END IF 
END 



121 



APPENDIX D. PRESSURE PROBE CALIBRATION CHART 



FT 

ill 



CALIBRATION DATA TOR TYRE DA & DAT 3-DIMENSIONAI. DIRECTIONAL TROBES 
t 7 ' IT] iiURTTTiiiTn: -TTnTnnanrj. Tn^n angle ^TnrrTi^iT^ETii^HinnHT:- 

-50* -40’ -30* -20* -10° 0* tio 1 + 20* + 30* +40* + S0*];ij 




wtiimowN 



PROBE TYrE:l>^-RS 
SERIAL NO.: ^>-2233-2 
DATE: k-3o-$& 

CALIB. BY: Vv tVe\r&- 
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APPENDIX E. REVX PROGRAM 



* OPEN A OLD FILE TO REVERSE THE X COORDINATE 
INTEGER A(300) 

REAL B(300),C(300),D(300),TEMP(500) 

REAL E(300)f(300),G(300),H(300),P(300),Q(300) 

OPEN( 1 JTLE='RESULTO.DAT') 
0PEN(2J 7 ILE=T1.DAT’,STATUS='NEW’) 

DO 40 K= 1,299 

READ(1 , 1 000) NO,X, Y,VEL, YAW, ALPHA J>T,CPTJ>S,CPS 
1000 FORMAT(I5,9F10.3) 

A(K)=NO 
B(K)=X 
C(K)=Y 
D(K)=VEL 
E(K)=YAW 
F(K)= ALPHA 
G(K)=PT 
H(K)=CPT 
P(K)=PS 
Q(K)=CPS 
40 CONTINUE 
CLOSE(l) 

N=299 
ICOUNT=N 
DO 10 1=1 ,N 
TEMP(I)=B(ICOUNT) 

ICOUNT=ICOUNT-l 
10 CONTINUE 
DO 20 1=1 JN 
B(I)=TEMP(I) 

20 CONTINUE 
DO 50 J=1,N 

WRITE(2,2000)J3(J),C(J),D(J)E(J)f(J),H(J),Q(J) 

2000 FORMAT(I5,7F10.3) 

50 CONTINUE 
CLOSE(2) 

END 
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APPENDIX F. XPLANE PROGRAM 



CHARACTER* 14 FNAME 

REAL Y(500),Z(500),YN(500),ZN(500),V(500),VXP(500), 

BETA(500),ALPHA(500),VX 1 (500),VY(500),VZ(500) 

C 

PI = 4.*ATAN(1.) 

WRITE(*,'(A\)')' INPUT DATA FILE? (D:FILENAME.EXT): ' 

READ(*,'(A14)’) FNAME 
OPEN(5JTLE=FNAME,STATUS='OLD') 

WRrrE(*,'(A\)')' OUTPUT FILE NAME? (DiFTLENAME.EXT): ’ 

READ(*,'(A14)’) FNAME 
OPEN(6,FILE=FNAME,STATUS=NEW , ) 

C 

READ(5,100) NDATA 
VXMAX = 0.0 
C 

DO 500 1=1 NDATA 

READ(5,101) I,Y(I),Z(I).V(I),BETA(I),ALPHA(I) 

C FIND CROSSPLANE VELOCITY' COMPONENT, VXP 

C CONVERT DEGREES TO RADIANS 

BETA(I) = BETA(I)*PI/180. 

ALPHA(I) = ALPHA(I)*PI/180. 

SB = SIN(BETA(I» 

SA = SIN(ALPHA(I)) 

VZ(I) = V(I)*SA 

VY(I) = SQRT(V(I)**2-VZ(I)**2)*SB 
VX1(I) = SQRT(VY(I)**2+VZ(I)**2) 

VXP(I) = VX1(I) 

C FIND THE LARGEST CROSSPLANE VELOCITY' COMPONENT 

IF(VXP(I).GT.VXMAX) VXMAX = VXP(I) 

500 CONTINUE 

NORMALIZE CROSSPLANE COMPONENTS TO A MAXIMUM VALUE OF 1 
DO 510 1=1 NDATA 
VXP(I) = VXP(I)/VXMAX 
0 CONTINUE 

— NORMALIZE CROSSPLANE VELOCITY TO A SPECIFIED FRACTION OF ONE 
— - GRID STEP FOR PLOTTING. ASSUME CONSTANT SIZE STEP. 

WRITE(*,'(A\)’)’ WHAT IS THE PROBE GRID STEP SIZE (IN.)?: ' 

READ(*,*) STEP 

WRITER, '(A)')' WHAT FRACTION OF GRID STEP IS MAX VELOCITY FOR’ 
WRnE(*,'(A\)')’ PLOTTING? (LIKE 0.5, ETC.): ’ 

READ(*,*) FRAC 
DO 520 1=1 NDATA 
VXP(I) = FRAC*STEP*VXP(I) 

520 CONTINUE 

C DIVIDE PLOTTING VELOCITY INTO COMPONENTS FOR COORDINATES 

DO 530 1=1 NDATA 

YN(I) = Y (I)+ V Y (I)/VX 1 (I)* VXP(I) 
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ZN© = za)+vz(i)^xiarvxpa) 

530 CONTINUE 

C STORE OUTPUT FOR PLOTTING BY "ARROW" 

WRITE(6,100)NDATA 

DO 540 I=1,NDATA 

WRITE(6,101)I,Y(I),Z(I),YN(I),ZN(I) 

540 CONTINUE 

100 FORMAT(I5) 

101 FORMAT(I5,5F10.3) 

CLOSE(5) 

CLOSE(6) 

STOP 

END 
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APPENDIX G. ARROW PROGRAM 



£*************** *********************************************** ******* 

c 

C THIS PROGRAM PLOTS DATA IN THE FORM OF VELOCITY VECTORS. 

C INPUT Y-Z POSITION, VELOCITIES, AND YAW AND PITCH ANGLES. 

C 

£**************************************************»****************** 
CHARACTER* 14 NAMEJTMAME 

C SF IS THE PHYSICAL SCALE FACTOR BETWEEN DATA AND PAPER. 

C 0.025 MM = ONE PLOTTER UNIT. 

WRITE(*,'(A)7 THE SCALE FACTOR SCALES THE PHYSICAL DIMENSIONS' 
WRITE(*,'(A)7 TO FIT THE PAGE.' 

WRITE(*,’(A\)')' DESIRED SCALE FACTOR? (FOR EXAMPLE, 0.5): ' 

READ(*,*) SF 

FACT = 25.4 * SF / 0.025 

WRITE (*,60) 

60 FORMAT (' ARROW HEAD LENGTH = ? (0.075 IS A TYPICAL VALUE) 1 ) 
READ (* ,*) HEAD 
HEAD = 25.4 * HEAD / 0.025 
WRITE (*,70) 

70 FORMAT (' ARROW HEAD WEDGE ANGLE = ? (30 DEG IS TYPICAL)’) 

PI = 4.0 * ATAN(l.O) 

READ (*,*) ANGLE 
ANGLE = ANGLE * PI / 180.0 

C "YREFZREF" IS THE INITIAL (REFERENCE) POINT RELATIVE TO WHICH 

C ALL DATA ARE PLOTTED. (0,0) IS ASSUMED. 

WRITE(*,'(A)')’ COORDINATES ARE ASSUMED TO BEGIN WITH 0,0 IN THE’ 
WRITE(*,’(A)7 FAR LOWER LEFT CORNER. IF THIS OK, TYPE 1. IF 
WRITE(*,’(A\)')' AN OFFSET IS DESIRED, TYPE 0: ’ 

READ(*,*) NCOORD 
IF(NCOORD.EQ. 1) GO TO 150 

WRITE(*,’(A\)7 DESIRED SHIFTED ZERO REFERENCE? (IE, -2., -2.): ’ 

READ (*,*) YREF.ZREF 
GO TO 151 

150 YREF = 0.0 
ZREF = 0.0 

151 CONTINUE 

C WRITE (*,200) 

C200 FORMAT (’ OUTPUT DATA FILE?: ’) 

C READ (*,*) NAME 
C CALL ZINIT (IPLOT JPORT J4AME) 

NAME = ’B:PLOT.DAT 
CALL ZINIT (1,1 .NAME) 

CALL ZVS(12.0) 

300 WRITE (*,’(A\)’)’ INPUT DATA FILE? (D:FILENAME.EXT): ’ 

READ (*,’(A14y) FNAME 

OPEN (5,FILE=FNAME,STATUS='OLD r ) 

C 

C THIS IS THE PROGRAM PROPER 
C 

READ(5,102) NDATA 
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DO 4001= l.NDATA 
READ (5,102) I, Y, Z, YN , ZN 
102 FORMAT(I5,4F10.0) 

C 

C THIS SECTION PLOTS THE TAIL OF THE ARROW (PEN UP) 

CALL ZPU 

IX = NINT((Y + YREF) * FACT) 

IY = NENT((Z + ZREF) * FACT) 

CALL ZPA (IX JY) 

C THIS SECTION PLOTS THE SHAFT OF THE ARROW (PEN DOWN) 
CALL ZPD 

IX = NTNT((YN + YREF) * FACT) 

IY = NTNT((ZN + ZREF) * FACT) 

CALL ZPA (IX JY) 

C THIS SECTION PLOTS THE HEAD OF THE ARROW 
CALL ZPD 

IF(ABS(YN-Y).GT.0.0001)GO TO 250 
ALPHA = -PI/2. 

IF(ZN.GT.Z) ALPHA = PI/2. 

GO TO 251 

250 ALPHA = ATAN( (ZN - Z)/(YN - Y) ) 

251 IF ( (YN-Y).GE.O.O )THEN 

IXU = IX - NTNT(HEAD * COS(ALPHA - ANGLE)) 

IYU = IY - NINT(HEAD * SEN(ALPHA - ANGLE)) 

ELSE 

IXU = IX + NLNT(HEAD * COS (ALPHA - ANGLE)) 

IYU = IY + NINT(HEAD * SIN( ALPHA - ANGLE)) 

END IF 

CALL ZPA (IXU.IYU) 

CALL ZPU 
CALL ZPA GX JY) 

CALL ZPD 

IF ( (YN-Y).GE.O.O )THEN 
IXU = IX - NINT(HEAD * COS(ALPHA + ANGLE)) 

IYU = IY - NINT(HEAD * SIN(ALPHA + ANGLE)) 

ELSE 

IXU = IX + NINT(HEAD * COS(ALPHA + ANGLE)) 

IYU = IY + NINT(HEAD * SIN( ALPHA + ANGLE)) 

END IF 

CALL ZPA GXU.IYU) 

C 

400 CONTINUE 

C CLOSE UP SHOP 

CALL ZPU 
CALL ZFIN1S 
CLOSE (5) 

STOP 

END 
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APPENDIX H. RUN MATRIX 







Test Run 


Grid #0 


Grid #1 


1 


Date 


Nov. 18, 1988 


Nov. 25, 1988 


Nov. 26, 1988 


2 


Start/End Time 


1526/2148 


0830/1821 


0143/1115 


3 


TwiA\vf 


62/66 


62/68 


62/78 


4 


p 0 (in Hg) 


30.14/30.19 


29.81/29.89 


29.97/30.06 


5 


PH 2 O (cm) 


7.2 


7.2 


10 


6 


Grid# 


0 


0 


1 


7 


Missile configuration 


B 


B 


B 


8 


AoA 


50° 


SOP 


50° 


9 


Measured dimension 


6.5"x8" 


3"x5.5" 


3"x5.5" 


10 


Step distance 


0.5" 


0.25" 


0.25" 


11 


Data file name 


R0B1B2O1.DAT 


R0B1A101.DAT 


R1B1A101.DAT 






to 


to 


to 






R0B1B214.DAT 


R0B1A113.DAT 


R1B1A113.DAT 


12 


Calibration file name 


CALDATT1 .DAT 


CAL0B.DAT 


CAL1B.DAT 






CALDATT2.DAT 


CAL0A.DAT 


CAL1A.DAT 
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